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I. INTRODUCTION 



The rapid increase of the structure function F2{x,Q^) at small x that is observed in 
DESY at HERA (see e.g. S ) has revived interest in the problem of the high-energy behavior 
of QCD amplitudes. In the leading logarithmic approximation it is governed by BFKL 
equation i i □ leading to a ~ behavior of F2{x) which is not far from the experimental 
curve. Unfortunately, there are theoretical problems with the BFKL answer which make 
it difficult, if not impossible, to use these leading logarithmic as a description of real high- 
energy processes. First and foremost, the BFKL answer violates unitarity and therefore it is 
at best some kind of preasymptotic behavior which can be reliable only at some intermediate 
energies. (The true high-energy asymptotics would correspond to the unitarization of the 
leading logarithmic results but this is a problem where nobody has succeeded in 20 years 
and not because of lack of effort.) 

Moreover, even at those moderately high energies where unitarization is not important, 
the BFKL results in QCD are not completely rigorous. Even if we start from the scattering of 
hard objects such as heavy quarks, then already in the leading logarithm approximation we 
obtain considerable contributions fruoj the region of small momenta (large distances) where 
perturbative QCD is not applicable □ □ . In other words, the hard pomeron which is believed 
to describe the observed small-x growth of structure function F2 interacts strongly with the 
soft "old" pomeron made from non-perturbative gluons. Therefore, it is highly desirable 
to have a method of separation of small- and large-distance contributions to high-energy 
amplitudes, and the starting point here must be a properly gauge-invariant formalism for 
the BFKL equation. 

In present paper we suggest a kind of gauge-invariant operator expansion for high-energy 
amplitudes. The relevant operators are gauge factors ordered along (almost) light-like lines 
stretching from minus to plus infinity. These "Wilson-line" gauge factors correspond to very 
fast quarks moving along the lines (see e.g. i ). It turns out that the small-x behavior of 
structure functions is governed by the evolution of these operators with respect to deviation 
of the Wilson lines from the light cone; this deviation serves as a kind of "renormalization 
point" for these operators. In this language the BFKL equation is simply the evolution 
equation for the Wilson-line operators with respect to the slope of the line. The gauge- 
invariant generalization of the BFKL equation turns out to be a nonlinear equation which 
contains more information than the usual BFKL equation — for example, it describes also 
the triple vertex of hard pomerons in QCD (cf. Q ). 

Asymptotic expansions (in large momentum limits) play a vital role in QCD. Cross 
sections (or amplitudes) in these limits simplify drastically, and one is thereby enabled to do 
calculations that would otherwise be impossible. The best established of these expansions is 
Wilson's operator product expansion for the T-product of two electromagnetic currents: 



Here the coefficients c„ contain all the singularities at x = 0, and the operators 0„ have 
no dependence on x. Taking the expectation value of Eq. (|l|) in a nucleon state and then 
Fourier transforming gives integer moments of the factorization theorem for deep-inelastic 
structure functions: 
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F2{xb, Q^) = Ci{xB, ® fi{xB, yU^ ".(/U^)) + . . . . 



(2) 



Here the parton densities fi{xB, /i^, cts(/^^)) are matrix elements of light-cone operators. The 
dots stand for the contributions of higher twist terms, i.e., terms damped by extra powers of 
l/Q^. xb is the Bjorken scaling variable xb = j'^V ■ <?, and /i is the renormalization scale. 

Both Wilson's operator product expansion and the factorization theorem can be expressed 
in terms of coefficient functions and operator matrix elements. This implies that a precise 
definition can be given to the quantities involved. In particular, there are contributions to the 
cross sections that come from the non-perturbative domain of large distances. The matrix 
element factors include these contributions, and their definitions include non-perturbative 
contributions. 

The renormalization scale /i has the qualitative effect of separating "hard" and "soft" 
contributions to the cross section. Integrals over soft momenta, those much less than /x, 
give suppressed contributions to the coefficient functions. Integrals over hard momenta, 
those much greater than yU, give suppressed contributions to the matrix elements. Roughly 
speaking the coefficient functions are given by integrals over large momenta Q'^ > p"^ > fi^, 
while the matrix elements are given by integrals over small momenta < /x^. The crucial 
property that enables calculations to be done easily is that the /i dependence is given by 
the renormalization group equations. One can set /i = 0{Q) in the coefficient functions. 
Both these and the kernel of the renormalization group equation can then be calculated 
perturbatively, in powers of as{Q)- 

Let us recall how the usual Wilson expansion helps us to find the dependence of the 
moments of structure functions of deep inelastic scattering. The essence is that instead of 
the dependence of the physical amplitude on (in the Euclidean region, which corresponds 
to the moments of structure functions), we study the dependence of matrix elements of local 
operators on the renormalization point /i. Consider the simplest Feynman diagram shown 



in Fig. |I|. 



q 



X 



0." 




FIG. 1. A typical diagram for deep inelastic scattering. 



At large q we can expand the current quark propagator in inverse powers of q 




(3) 
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where n*^ term of the expansion corresponds to the n*^ moment of the structure function. 
Unfortunately, after the expansion the loop integrals over k become UV divergent so we must 
modify our Taylor series (|^) somehow. To this end, we note that each term on the right-hand 
side of Eq. corresponds to the matrix element of a certain quark operator of the type 
i/j[dfj_)"''jiy'ip] the UV divergence reflects merely the large dimensions of these operators. It is 
well known how to regularize these UV divergences — we must introduce the regularized 
operators normalized at some point fi and expand our physical amplitudes in a series of 
these regularized operators. Roughly speaking, each term on the right of Eq. (P) will be 
integrated over k only up to /c = /i. The dependence on /i will be cancelled : in the next 
order in the coefficients of the Taylor expansion will be modified also — they will contain 
terms ~ In ^ which will cancel the dependence of matrix elements of the (renormalized) 
local operators on /x^. In the leading logarithmic approximation we can simply take /i^ = 
and the dependence of moments of structure functions on Q'^ will refiect the dependence of 
the matrix elements of the operators ^V^^...V^„_i7^„V the normalization point. This 
dependence is given by renormalization-group equation (see e.g. the book 13 ). 

Summarizing, in order to find the dependence of the structure functions of deep inelastic 
scattering at large we perform the following steps: (i) formally expand in inverse powers 
of Q^, (ii) regularize the obtained UV divergent matrix elements of local operators, and 
(iii) write down (and solve) the evolution equation with respect to the normalization point 
/i. In the original Feynman diagrams for structure functions of deep inelastic scattering 
the photon virtuality plays the role of a "physical" cutoff for the integrals over loop 
momenta. After expansion in powers of these integrals became UV divergent; by 

adding counterterms in the usual way we introduce an "artificial" cutoff /x^ for these loop 
integrals. Now, instead of studying of the behavior of the original Feynman diagrams, we 
should trace the dependence of the matrix elements of the operators on the cutoff /x^. This 
is a lot easier thing to do because it is governed by renormalization group. 

Now, we would like to generalize these ideas for high-energy scattering. In order to find 
the high-energy behavior of a certain physical amplitude (say, the structure function of deep 
inelastic scattering at very small x), we will perform the same three steps: (i) formally 
expand the amplitude at large energy s — after that we will have the divergences in the 
longitudinal integrals, (ii) regularize these longitudinal divergences by introducing a certain 
cutoff, and (iii) find (and hopefully solve) the evolution equations with respect to this cutoff. 
As in the case of Wilson expansion, the dependence of the relevant matrix elements on the 
cutoff determines the high-energy behavior of the original amplitude. In the subsequent 
three Sections we will perform these steps. In the Appendices we present the shock-wave 
picture of high-energy scattering in the virtual photon frame. 



II. HIGH ENERGY LIMIT 

As an example, let us consider the high-energy behavior of the forward scattering ampli- 
tude for virtual photons in the region where s = {pa + PbY ^ P%P%- 
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AipA,PB) = -^eX^!^v I d'xd'ydhe''P--^^'P--y{0\T{f{x + z)f{z)j^iy)f{0)}\0). (4) 

(Only the connected part of the Green function is used, and the vectors and are the 
polarizations of the photons.) For simplicity, we assume that the virtualities of photons are 
negative, since then our amphtude (H) will have only one discontinuity corresponding to the 
total cross section of virtual photon scattering. (At s ^ ^ this will be the cross 
section of deep inelastic scattering from a virtual photon at small x.) A typical graph is 
shown in Fig. |^. Our aim is to obtain the leading contribution (in powers of s) from graphs 
for the amplitude, and it is well known that the leading large-energy asymptotic behavior (s 
times logarithms of s) corresponds to diagrams with gluon exchanges. 




FIG. 2. High-energy scattering of virtual photons. 

To orient the reader in our subsequent technical treatment, we first explain the qualitative 
features of the process in coordinate space. Suppose that we view the process as one of the 
incoming photons, A say, travelling through the color field due to the other. Because of 
time dilation and Lorentz contraction, each beam particle may be viewed as a collection 
of fast-moving, point-like objects distributed over the transverse plane. The probability of 
a large momentum transfer Q is of order l/Q"^, so that the partons are to be regarded as 
travelling along straight lines while they are crossing the non-trivial part of the field. This 
is, of course, the parton model. In a lowest order approximation, such as Fig. 0, the partons 
in question are given by a single quark-antiquark pair. The photon has fluctuated to a state 
of a quark-antiquark pair, and this state is almost unchanged while the pair traverses the 
field, since the time-scale for the evolution of the state is much longer than the time to cross 
the field. 

Another view of the same situation can be obtained in the rest frame of the beam. The 
field of the other particle is Lorentz contracted, time dilated (and intensified). It looks like 
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a shock wave of width L ~ \pW/s, and the quarks of the beam cannot make a significant 
movement in the transverse direction during that time. Here, L is the typical length scale 
of the field. 

From either viewpoint, we see that the situation is one in which some version of the 
eikonal approximation is valid. That is, the effect of the field on the state of the quark- 
antiquark pair is given by integrals over the gluon field along the (straight-line) trajectories 
of the partons. Indeed, we will see that what we need are exactly path-ordered exponentials 
of the gluon field. Unlike the simple eikonal approximation, the field affects both the phase 
and the color orientation of the state. 



A. Amplitude as integral over gluon field 

As usual, at high energies it is convenient to use a decomposition in Sudakov variables. 
For the momenta we write the standard formula 

= + (5) 

2 2 

where ~ — and P2 — Pb ~ ^Pa light-like vectors close to pa and pb 

respectively. (Then gf^^, = |(pi^p2i/ +PluP2^l) + gj[y) These variables are essentially identical 
to light-front coordinates, a = p+/^/s, (3 = p^j^. For the coordinates we use a slightly 
different form 

= lz,pi^ + lz,pi^ + z1, (6) 

s s 

where 2;. = z^pj* and = z^p'^'^ One advantage of these coordinates is their simple scaling 
properties when we take the high energy limit, as in Eq. (|T^, below. The factors 2/s in 
the formula for the components, Eq. (^, avoid extra factors of s in the combination p ■ z 
=apZ, + PpZ^ -p±- z±. 

The Jacobian of the transition to Sudakov variables is s/2 so that 

J d'^z = - J dz,dz^:(fzj_ , J d^p = ^ J dapd(3pd'^p±. (8) 

To put the scattering amplitude (^ in a form symmetric with respect the top and bottom 
photons, we make a shift of the coordinates in the currents by {z,,0,0±) and then reverse 
the sign of z,. This gives: 



" Sometimes, however, we shall use also "covariant" coordinates - Sudakov variables: 

= up'^ + vpt^ + z'l, (7) 
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A{pa,Pb) = -le^e^efe^'^ Jd^Adz, Jd^xd^ye'P^-''-^'P^-y 

(0|T{j'^(x., x, + z,,x^ + z^)fiO, z,, z^)f{y, + z„ y,, y±)f{z., 0, 0±)}\0). (9) 

(We remind the reader that only the connected part of this Green function is taken.) 

It is convenient to start with the upper part of the diagram, i.e., to study how fast quarks 
move in an external gluonic field. After that, functional integration over the gluon fields will 
reproduce us the Feynman diagrams of the type of Fig. ||: 

Apa,Pb) = -ie^^etefef-jd\^M-^ j VAe'^'^^MeiiiV) 

[iJdz.l d^x e'P^-^' {Tf{x„ X, + z,, x^ + zx)j"(0, z,, z^)) ^ 

{-Jdz.J d'y e^--y {Tj^iy. + z., y., y±)f{z., 0, 0^))^} , (10) 

where 

{Tj,{x)j,{y))A = jVtfjVtpeSi^'^) at A = 0, (11) 

and S{A) and S{ip, A) are the gluon and quark-gluon parts of the QCD action respectively. 
det(zV) is the determinant of Dirac operator in the external gluon field; it gives the effect of 
quark loops in Fig. |[ The subtraction in Eq. ([ll|) removes the disconnected graph. 

The arrangement of the integrals in Eq. (|l^) arises from a choice to construct amplitudes 
that have all momentum conservation delta functions removed. The integrals over x and 
y set the momenta of the outgoing photons to be pa and pb- The integral over z^, sets 
the P component of the incoming photon momentum on the top bubble to be equal to the 
corresponding component of the outgoing momentum, f3p^ ~ p\/ s. The a component of the 
incoming photon to the top bubble is the corresponding component for the outgoing photon 
minus whatever a component of momentum the external field happens to provide. A similar 
statement applies to the z, integral. The z_i_ integral enforces zero transverse momentum 
transfer at one end, and we leave it as the outermost integral in order to emphasize that 
we wish to treat transverse momenta symmetrically between the upper and lower quark 
loops. There remains the functional integral over the gluon field, after which momentum is 
conserved. Before this is performed, there is no conservation of momentum, since the gluon 
field is position dependent. 



B. Fast moving photon in external gluon field 

From Eq. (0), it is clear that the amplitude for the upper part of the diagram in Fig. |^, 
describing a virtual photon with momentum pa flying through the external gluon field A^, 
is given by the following expression (see Fig. ^): 
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FIG. 3. High-energy scattering of virtual photon from external field. 



Jdz^ J d^xe'^^'''{T{jf,{x,, + z*, x± + z±)>(0, z^, z±)})a 

2, r i 

= - ^ef / dz^d^xe'^^'"' Tr 7^([x,, x^ + z^,x± + z± — 0, z^, z± 
s ■ J f 



7,,((0, z^, zx. ^ x„ x^ + z^,x± + z±Jj 



—Same at A = 0, 



(12) 



where we have used Schwinger's notations for the propagators in an external field, and Cj 
are the quark charges. 

In Schwinger's notations we write down formally the quark propagator in the external 
gluon field A.,( clS db matrix element of the inverse Dirac operator 



G{x,y) 



X 



where 



d 



X 



(13) 



{{x\y)) = 5^'\x - y), {{x\p,\y)) = - y), ((x\AM = A,{x)S^'\x - y). 

(14) 

Here are the eigenstates of the coordinate operator A'lx)) = (normalized according 
to the second hne in the above equation). From Eq. (|14D it is also easy to see that the 
eigenstates of the free momentum operator p are the plane waves = / d^x e~'^P'^\x)). 
It should be clear from the context whether the vectors are eigenstates of momentum or 
position. Note that the states or are functions of /our- vectors {x^ or p^), unlike the 
actual quantum mechanical state vectors of the field theory. 

Thus, for example, the first term of the expansion of the propagator (^) in powers of 
external field is the free propagator 



16n^ f 



(15) 
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We are treating the gluon field as a matrix in the fundamental representation of SU(3): 

A,{x) = J2A",t^- (16) 

a 

Then the quark propagator, Eq. (|T^), is a matrix in both color and spinor space; the = 'j^p^ 
is implicitly multiplied by a unit color matrix. 

Now let us Fourier transform Eq. (|12|) over z±, so that it is a function of q± instead of 
z±. Going to Sudakov variables (S), we have: 



s2 2 /• d^k d^p 
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Tr {74(/^|^|fc-p))7.((fc-PA-p'|^|A:-p4} 

-Same at A = 0. (17) 

Here, {ah)±_ denotes a scalar product of transverse components of vectors a and h. In this 
expression, the quark- loop momentum is = auPi + [3kP2 + while = apPi + PpP2+Pi_ 
and p''^ = Up'Pi — (3pP2 — P± + ci± are the momenta entering the two quark lines from the 
external field. Notice that the quark propagators do not conserve momentum. We have 
enforced conservation of the f3 and the transverse components of momentum by our choice 
of external momentum, while conservation of the a component of the momenta will only be 
true after the functional integral over the external gluon field to form the complete amplitude, 
as in Eq. (^). 



C. Regge limit 

Now, we must formally take the limit s — > oo in this expression. We will do this for 
a fixed external field. The Regge limit s ^ oo with p\ and p% fixed corresponds to the 
following rescaling of the virtual photon momentum: 

PA = \P?^ + —^—P2, (18) 

2Xp\ ■ p2 

with pb fixed. This is equivalent to 

p^ = Xpf\ p^=p^^\ (19) 

where p^i^ and p^2^ are fixed light-like vectors so that A is a large parameter associated with 
the center-of-mass energy (s = 2\p^i^ ■^2°^)-^ 

''The method we are using is a version of the methods used in iH . 
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Next, let us look at how the Sudakov variables in Eq. p^ ) scale with A. In general 
when treating an asymptotic limit of some Feynman graphs, there will be a number of 
different regions of loop-momentum space that contribute. It is quite a complicated problem 
to disentangle these. However, we have chosen to start with the asymptotics for a fixed 
external field. So initially we do not have to conern ourselves with the problem of multiple 
regions. That problem arises at a later stage of the argument when we perform the the 
functional integral over the gauge field. 

The limit we are taking is A ^ cxd with the gauge field A fixed. 

First, we shall see below from the explicit form of integral that the important values of 
the variables for the quark loop attached to photon A satisfy ~ 1 and ~ ^. Such 
momenta are obtained by boosting from A = 1. With this scaling, scalar products of quark 
momenta, and the measure d'^k are independent of A; this is a consquence of boost invariance. 

Moreover, the important values of momenta transfered from the gluon field obey ap ~ 
j,Pp ~ 1, since A{ap,Pp,p±) = A^'^^Xap, Pp,p±) 

{a, f3, p^) = f d^x A{x) e^°pr-^+iM°*-^-i(p^)i (20) 



is a function independent of A f (Hereafter the A{p) denotes the Fourier transform of the 
field A{x)) So, we must compute the behavior of the quark loop (|1^ in the region 

ak ~ ~ ~ 1 {-Pa) 

«P ~ ^, /5p ~ 1, pi ~ 1 i-pl). (21) 



(We shall see that k\ ~ p\ from the explicit integral (|36D and that the characteristic pi of 



the external field is determined by the characteristic scale of the source of this field, which 
is the virtuality of the target photon p^.) 



D. Quark propagator in external field 

As a first step, we will find the quark propagator in the external field (see Fig. ^). In 
the limit we are considering, we must recall the well known fact that at high energy we can 



■^The scaling for A applies before the functional integral over A. After the integration over A, we will get 
contributions from ap ^ 1 and from ap — s- oo. The first region corresponds to higher-order corrections to 
the quark loop; these are just like higher order corrections to the Wilson expansion. The second region 
corresponds to UV divergences in these same higher-order corrections. In both cases subtractions must be 
applied; we treat this as a separate issue. 
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replace g^^, for the gluon propagators connecting quark lines with very different rapidities by 
^Pi^P2u- Thus, we can change the factors 7^A^ for the interaction to -J)^A,^ correct to the 
leading power of A (or s). This gives: 



-9 



^9 



s + ie 
,2 f d'p' 



{k — pY + ie 



s J 167r4 /c2 + 



le 



{k — p')2 + ie 



{k — pY + ie 



+ ..., (22) 
where the dots stand for further terms in the expansion in powers of the external field. 




FIG. 4. Quark propagator in the external field. 

Let us start with the first nontrivial term ~ A, shown in Fig. ^3. From the eq. ( ppj ) it 
follows that 

A(ap,/3p,p±) = XAf\\ap,f5p,pA_) (23) 

(here Ao = A^p^^). Now it is easy to see that in the limit A ^ oo the Fourier transform of 
the external field A{p) is proportional to 6{ap) (we assume that the Fourier transform of the 
external field A^^\p) eq. ( pOD decreases at infinity). The coefficient in front of the (5-function 
can be figured out from the following formula 

tg J ^A.{ap,l3p,p±) =tg Jd\±e-'^'^^^ j du A,{upi + x,_)e^'^-^ . (24) 
so the first two terms of the expansion of the propagator (^) reduce to 
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1 

f 



k — pj 



-g- 



-iip,x)j 



du A,{upi + x_L)e2 



{k — pY + 



(25) 



As we will see below, the integrals in Eq. (p!7D will force j3p to be of order 1/A, so that we 
should set /5p = in Eq. (p^. Eq. (p^ ) is the first order of the expansion of a path ordered 
exponential whose precise definition is given in Eqs. ( p9D and (|31| ) below. 

Next, consider the third term on the right-hand side of Eq. (^) shown in Fig. ^c. In the 
region (^), the second propagator of this term (times the factors on each side) reduces 
to an eikonal denominator: 



i>2 



{ak - a'p)^^ + (^ - 



(ftfc - a'){Pk - P')s -{k- p')l + ie 



(afc - ap(/?fc - p'p)s -{k- p')l + ie 

^2 o, ^ • • (26) 

-13'^ + leak ^ ' 



Furthermore, we can neglect compared to in the free quark propagators in Eq. 

f|22|). 



Again, since characteristic — at A ^ oo this contribution will be proportional 
to S{ap). In order to find the coefficient in front of this S function we shall integrate the 
propagator over ap. As we will see below, when we do the integral over f3k, both the quark and 
antiquark lines are restricted to being forward moving, i.e., < < 1. Then the integrals 
of the gluon fields in this third term reduce to the second order term in the path-ordered 
gauge factor: 



dap da'p d/S'p d^p'^ ~ 



2ti 27r 27r 47r2 



A,{a' [3' p'^ 



1 



-A,{ap- ap,[3p- (3'p,p^ -p^] 
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d\^ e^'^P''^^' 



-g'^ / d'^x±e 



-P'p + le 

dv j du<d{v - u)A,{vpi + x±)A,{upi + x±) e^''^''" 
dvfdue(.-u)A.(.p. + .M.(up. + -.)- (27) 



-iip,x)j 



In the last line, we have again used the result, to be demonstrated later, that Pp is of order 
1/A. So, the expansion (p^ takes the form: 



f 

IGtt^ 



5^'\p) 



g-2TT6{a. 



J2 



J d'^x± e'^^P'""^^ J du A,{upi + x±)e^ 



sf^pU 



{k — pY + 
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J (fx± e~^^P'^)i- j dv j duQ{v - u)A,{vpi + XA_)A,{upi + 



+.... (2^ 



{k — + ie 



We now express these and all the higher terms of the expansion of the right-hand-side of 
Eq. (|22|) in terms of path ordered exponentials of the gluon field. Let us use [x, y\ to denote 
the path-ordered gauge factor along the straight line connecting the points x and y: 



\x, 



_ p^ig du{x-yY A^{ux+{l-u)y) ^29) 



Then it can be demonstrated fairly easily that further terms of the expansion of the right- 
hand-side of Eq. (|2^ ) in powers of A, will reproduce the subsequent terms in the expansion 
(in powers of A,) of the following gauge factors:^^ 

[U ~l]{p^) = U{p^)-A7,'5^'\p^), 

P - = U\p^) - 47r25(2)(p^), (30) 

where U{p±) and U'^ipx) are Fourier transforms of gauge factors along lines extending to 
infinity in both directions: 

U{xi_) = [oopi + x±, -oopi + 

U''{x±) = [-oopi + x±,oopi + x±]. (31) 

These correspond to quarks moving across the external field with the speed of light. The 
relevance of these gauge 

Thus we finally have the propagator of a fast-moving quark with < < 1: 



This is valid to the leading power of A, when A is large and the external gluon field A is 
fixed. A similar formula is valid for the antiquark propagator. 



'^The form of our notatfons, [[7 — 1], [W — 1], reflects the fact that at /3p = 0, these gauge factors are simple 
path-ordered exponentials along an infinite line, but that the zeroth term in the expansion in powers of 
gauge field is missing. 
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E. Impact factor 



Let us rewrite the expression (|T^ for the upper part of the diagram in Fig. ^ using the 
above formula for quark propagator (|32D, and the corresponding formula for the antiquark 
propagator. We obtain 

-i:«?/^e(i>"^>o) 



Tr <i7A. 



k2 + ,e [(A; -PA + q)l + ^e] [(A;p^)2 + le] ^ ^ 

[A;i + ^e][(A; - qf + ze] {k - p^)^ + ^^"^ '^^^^^^ 

[(A;-p^-p + g)2 + 2e][(A;-p^)2 + ze] J^^^ J' 



(33) 



where we used the notation ak = 1 — «fc- Now we can use contour integration to perform 
the integrals over and f3p. It is easy to verify that the dominant contribution arises when 
both these variables are of order a squared transverse momentum divided by s, i.e., of order 
1/A. It is easy to see that the linear terms in U and f/^ cancel in Eq. (^) so after some 
algebra one obtains the final answer in the following form*^ 

= T.^'J ^I,Ap±,q±)TT{U{p^)U\q^~p^)}, (35) 



'^A more careful analysis performed in Appendix |^ shows that the Wilson lines U and are connected by 
gauge factors at infinity so 

TTU{x^)UHy±) (34) 
lim Tr[-Api +x±,Xpi +x±][Xpi +x±,Xpi +y±][Xpi +y±,-Xpi +y±][-Xpi + y±,-Xpi +x±]. 

A — >oo 

The gauge factors connecting the end points of the eikonals U and reduce at infinity the gauge factors 
made from pure gauge fields so the precise form of the contour connecting the end points of Wilson lines 
does not matter. 
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where I^u{p,q) is the so-called "impact factor": 







s[(/c + + (gi J' 

Contrary to appearances, the impact factor is independent of s (and hence of A). The easiest 
way to see this is to observe that by a boost of the coordinates used in Eq. (|19|) we may 
obtain the large s limit by scaling But in Eq. (|3^) p2 only occurs in the combination 
V2IS. 

When the photon indices \i and v are transverse, we obtain the following explicit expres- 
sion for l^y. 

T I \ I da da' r „2 / / 2 / 2\ -~|-^ 
lMP±,q±> = -- — —\P^aa' + {q^a-p^)aa\ 
I Jo Zn Jo 2,71 ^ J 

|(1 — 2aa — 2a' a' + 8aaa'a')P]_g^,y + d>aaa'a' P±^P^y 

+2g^yq]_aa{l - 2a') - 4aa(l - 2a)Q;'Pj_^g_Li/|, (37) 

where P± = p± — q}_a. At g_L = this result agrees with . 

Note that the limit A — * 00 enforces the vanishing of the total (3 argument of the gauge 
factors U (and U'^) in Eq. (|35|), whereas the individual A, fields forming this gauge factors 
may have nonvanishing (3's. It is instructive to write down the final formula for the quark 
propagator in this case: 



f 

— ^K)tIt^ {P - l](p±)e(a,) - [t/t - l]{p^)Q{-a,)) . 
s k]_ + le ^ ^ [k — pY^ + le 



It is worth noting that the form of the answer (^) — free propagator ® eikonal factor ® 
free propagator — is due to the shock-wave structure of the external field at large energies 
(see Appendix ^ . 

Let us also present the result ( P5| ) in the transverse coordinate representation. One has 
for the forward scattering: 



--T.^iJ d^^i. j d?z^ itui^A.) Tr{f/(xx + z^)U\z^)], (39) 
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where the impact factor in coordinate representation has the form: 



1 dada' 



Ana'a' \ 



p\aa 
a'd' 



-2g^^{l - 2aa){l - 2a'a)Ki 



aa 



+ 1 



aa 



g^u{l — 2aa — 2a a + 9>aaaa) + 9>aaaa 



— I —I ^IJi^V 



aa 



(40) 



where if^^ipi.) = -^/fjy(p±,0) and Kn{z) is the McDonald function. 

Formula (|39|) describes a quark and antiquark moving fast through an external gluon 
field. After integrating over gluon fields (in the functional integral) we obtain the virtual 
photon scattering amplitude ([TT|) . It is convenient to rewrite it in the factorized form: 

'^'P'--I^{p^){{TT{U{p^)U^{-p^m. (41) 



Apa,Pb) 



E 



47r2 



where I'^{p±) = e^e^I^uiP±)- The gluon fields in U and W have been promoted to operators, 
a fact we signal by replacing U by U, etc. The reduced matrix elements of the operator 
Tt{U{p±)U^ {—p±)} between the "virtual photon states" are defined as follows: 

{{Tt{U{p^)U\-P^)})) = /dVe-^(^^)-((Tr{?7(x^)f/t(0)})) 



{{Tt{U{x^)U^{x'^)})) = - d'z6{z, 



{0\T{TT{U{x^)U\x'^)}j^{y + z)f{z)}\0). 



(42) 



It is worth noting that for a real photon our definition of the reduced matrix element can 
be rewritten as 



,PB\TT{Uix^)U\x'^)}\e',pB + f3pB) = 2nSi(3){{TT{U{x^)U\x'^)})), 



(43) 



where e and e' represent the polarizations of the photon states. The factor 2tt6{i3) reflects the 
fact that the forward matrix element of the operator U{x±)W{x'j_) contains an unrestricted 
integration along pi. Taking the integral over jS one easily reobtains Eq. (0). 

Our expression Eq. (|35D represents the upper part of the graph as a numerical factor 
times a function of the gluon field. The result is independent of what we chose to put in as 
the lower part of the Green function. Thus we may say that this formula is correct in the 
operator sense: 



d^x / dS5(;z.)e-*(«'")^e*P-4-^T{j^(x + z)j^{z)} 



E^ 

flavors 



^I>±,<l±)^r{U{p)Wiq-p)}, 



(44) 
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where the operators U and ?7"l" are given by the same formulas (|3TD with the substitution of 
the external field A by the field operator A. (We continue to use the ( ' ) notation for the 
operators in order to distinguish them from the corresponding expressions constructed from 
external fields.) 

This formula is a bit misleading, since the derivation assumes that the gluon field only 
has Fourier components that obey ftp <^ 1 and (3p < 1. However, we expect that Fourier 
components that do not obey this condition, in particular ~ 1, will effectively give higher 
order corrections to the coefficient J^. 

This is the first term in an expansion in powers of A at large A. In Eq. (44), J^^, has the 
same status as a Wilson coefficient: it is a numerical coefficient that multiplies an operator. 
However, unlike the case of the Wilson expansion, the coefficient is not a pure ultra-violet 
quantity; we plan to express it as yet another operator matrix element. 

Unfortunately, matrix elements of the operators TT{U{p)W{q — p)} ordered along the 
light-like line will have a longitudinal divergence in Feynman integrals. This is rather like 
the Wilson OPE where the matrix elements of the (unrenormalized) local operators will 
have UV divergence in the integrals over loop virtualities. In the next Section we will 
introduce "regularized" eikonal-line operators U and W which will be the analogs of the 
local renormalized operators for high-energy amplitudes. 



III. REGULARIZED WILSON-LINE OPERATORS 

In the previous Section we have found that the formal high-energy limit of the virtual 
photon scattering amplitude is described by a matrix element of a Wilson-line operator 
( Pip ordered along a light-like line. However, a matrix elements of such an operator has a 
longitudinal divergence. 

We will now explain how the divergence arises and how to treat it, with the aid of a low 
order example. 



A. General structure; divergences, subtractions 

Originally we had graphs of the form of Fig. |^. For the present part of our argument, 
let us choose all transverse momenta to be of some given order of magnitude. Call this 
magnitude m. (Finally we will see that all the integrals over tranverse momenta converge 
on scales of order of photon virtualities). Then the operator factorization Eq. ( PD applies 
as it stands when all the gluons have a -C 1. The longitudinal momenta are then restricted 
to give unsuppressed contributions only when they are not too big: < m^. (This 

last statement is actually in need of proof.) We are using a Sudakov representation for the 
momenta — Eq. (^). 

When we consider the integral over all the longitudinal gluon momenta, we can partition 
the graph into factors ordered from top to bottom. The a's are strongly ordered between 
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the different factors, with the largest values at the top. (We will not present the proof that 
configurations with strong reverse ordering between two factors are power-law suppressed.) 




FIG. 5. First correction to a lowest order graph for high-energy scattering. 

Let us add one extra gluon rung to a lowest order graph, so we have Fig. |^, and let us 
write the graph as 



G = J d'plip), 



(45) 



where is the momentum flowing from one of the gluon lines into the upper quark loop. The 
procedure we explained in the previous section gives us the asymptotics for the integrand 
when m^/s < <^ 1, where, as we defined earlier, m represents the typical scale of 
transverse momenta. So the factorized form is in an integral 



d^p Asy I{p), 



(46) 



with Asy / being of the form of the impact factor times the integrand for a graph for 
((Tr{f7(px)f/t(-p^)})). 

Then we write the graph. Fig. ^ as 



G= fd^p Asy I{p) + f d^p 



I{p) - Asy I{p) 



(47) 



The first term is the lowest order impact factor times a correction to the operator. The 
second term has its -C 1 behavior subtracted off, and so the dominant contribution to 
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the integral in this term is from of order unity, or bigger. This term should be treated as 
giving a higher order correction to the impact factor, as we now explain. 

Suppose that we have proved the factorized formula, Eq. (|41|), in general. Consider its 
expansion in powers of the coupling. The first term on the right of Eq. (^) is a contribution 

to 

lowest order impact factor x next order matrix element, (48) 
while the second term is a contribution to 

next order impact factor x lowest order matrix element. (49) 

But, as we will see shortly, the integral over Asy„^^^ I{p) has a divergence as ap oo, 
since the replacement of I{p) by Asy /(p) removes a convergence factor provided by the quark 
loop. (The approximations used to derive Eq. (^TJ) are only valid when ap <^ 1.) We must 
therefore redefine the operator U{p±)W{—p±) so that it has no divergence. Ideally we would 
like to do this by some kind of generalized renormalization procedure. But for our discussion 
we will find it sufficient to change the line along which the path ordered exponential is taken. 

The structure of Eq. (^^ and the arguments that we will need are completely analogous 



to those for the ordinary operator product expansion. However, it is important to realize 
that the divergence we are concerned with is not a conventional ultra-violet divergence. Thus 
the methods used for the operator product expansion need to be generalized. (The operator 
indeed has ultra-violet divergences, in certain graphs. These are associated with p± ^ oo 
behavior, and constitute a relatively trivial problem.) 

The decomposition of the amplitude into the impact factor times matrix element has 
a very illuminating (although qualitative) interpretation in terms of functional intergral 
representation for the amplitude. It corresponds to the decomposition of the functional 
integral into a product of two integrals - over the (quark and gluon) fields with large light- 
cone fraction a ~ 1 and over the fields with small a ~ ^ (which corresponds to fields that are 
not scaled with A). More presicely, we choose a such as a -C 1, g'^lna <^ 1 {a is independent 
of A) and separate the functional integration over the fields with lighcone fraction a either 
greater or lesser than a. First , we perform the integration over the a > a fields which scale 
with A and it yields impact factors times the Wilson-line gauge factors constructed from 
"external" fields with small a < a and on the second step the remaining integral over these 
small-a fields will give us the matrix elements of the Wilson-line operators. (In the leading 
logarithmic approximation these Wilson-line operators still correspond to the slope || pa 
since we make no difference between In ^ and In^). So, the impact factor is given by a 
functional integral over the a ~ 1 fields in the external small-a fields which technically is a 
series of diagrams in the external field. The leading-order impact factor calculated in Sect. 2 
is the simplest of such diagrams. In the next order in coupling constant we will have more 
complicated diagrams with large-a gluon fields as shown in Fig. 5. Unfortunately, there is 
no consistent quantitative decomposition of the functional ontegral into product of a > a 
and a < a integrals which goes beyond leading logarithmic approximation. So, at this point 
we are forced to return to the original logic of the operator expansion and define the impact 
factor as the coefficient function in front of the (Wilson-line) operator by comparing the 
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matrix elements of the T-product of two currents and of the Wilson-hne operators. If we 
knew that the expansion goes in terms of Wilson lines beforehand and our purpose was just 
to calculate the coeficients, it would be enough to compare these matrix elements between 
two (or four) real gluons. But since we want to prove that the gluon operators assemble in 
Wilson lines we must compare these matrix elements between an arbitrary number of real 
gluons - i.e., in external gluon field. So, again the impact factors are given by the diagrams 
in the external gluon field but the interpretation now is different - the external gluon field is 
a convenient way to represent many-gluon states between which we must take the operator 
expansion in order to determine the coefficient functions (impact factors). Maybe if the 
correct gauge-invariant way to separate functional integrations over large and small distances 
will appear some day it would possibly make the two interpretations of the same diagrams 
in external fields equivalent. 



B. Loop corrections to Wilson line operators 

Consider the example of the one-rung ladder diagram shown in Fig. p. 




FIG. 6. Typical diagram for the matrix element of the operator Tr{f/C/t} between "vir- 
tual-photon" states. 

The corresponding contribution to the matrix element ((Tr{f/(p_|_)f/"''(— p_|_)})) has the 
form: 



dap da'p d/S'^ d'^p'j_ 
'2^ J 27r 27r 2tt Att'^ 



piia'pP'pS -p'l + iey [-{ap - ap/^^s - {p - p')l + ie] 
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up to the trivial color factor Nc{N^ — 1). Here the momenta are defined in Fig. 

r 

tex, and 



ij,ua{p, -p',p' - p) = -{p + P')a9t,u + (2p' - p)f,gua + (2p - p')u9ai^ IS a three-gluon ver- 



1 



£ /If 



1 



1 



-7" 



(51) 



is the quark loop shown in Fig. ^ In this section we also omit for brevity the trivial factors 




FIG. 7. Quark bulb. 

due to the electric charges of the quarks. 

In writing Eq. (SD), we have assumed that rapidity of the gluon rung is much larger 
than that of the quark loop. This accounts for the indices on the 3-gluon vertices. We used 
Feynman gauge and substituted (7^^ by ^Pi^P2ri which is valid for gluons connecting lines with 
very different rapidities. 



C. Calculation of divergences 

It is easy to see that the integral over a in Eq. (|50D is logarithmically divergent. At first 
sight, the divergence appears to be linear, since 

lr,.'^(p, -p')r,.Ap, -p') = {a, - 2a'^)l3'^s - {p + (52) 

but careful analysis carried below shows that /^^ is ~ ^ at large ap.) The only gauge- 
invariant way to regularize this divergence which we have found is to change slightly the 
slope of the supporting line (as it was done in ref.0 for the case of Sudakov formfactor). 
We define: 

U'^{x^_) = [oop'' + x±, —oop'' + x±] , W'^{x±) = [—oop^ + x±, oop^ + x±] (53) 

where 
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Pi + CP2 



(54) 



so that at C <^ 1 the operators ([53| ) are ordered along a shghtly non-hght-hke hne. 

Now let us demonstrate on our example that changing of the slope of the line according 
to Eq. ( ^41) does regularize the longitudinal divergence in matrix elements of the opera- 
tors U. It is easy to see that the changing of the slope of the line according to Eq. (|5^) 
leads to the substitution p = api + p± — > p = api — CttP2 + P± in the diagram in Fig. ||. 
Therefore, we obtain the contribution of this diagram to the matrix element of the operator 
{{Tt U'' {p±)W'' {q± — p±))) in the following form: 



'2^ 



dap d p'j_ 
27r 167r4 



[(ap-2a')(3's-ip + p')l]^^ip') 



(Cap +pI- ief {a'^ii'pS - p'\ + itf \-{ap - a')(«pC + i^'p)s - (p - p')l + «e] 



(55) 



As we shall see below, the logarithmic contribution comes from the region y ^ 3> ctp 3> 

a'p 1 ^ jS'p ^ Pp = —(ap ~ V In this region one can perform the integration 

over (3'p by taking the residue at the pole — («p — a'){ap( + (3'p)s — (p — p')± + it , and the 
result is ^ 

s 



r dap 


da'p 


f d'^p'i. r 


1 27T 


271 J 


47r2 y 



6(ap > Op > + 6(0 > a'p> ap) 



{pl+p'l-alCs/2) ^''yappi 



Pa 



\ap - a'p\ (Cap + pi - ief [^(p 



P^Tl^pI 



(56) 



Here we have used the approximation that ^ a^: The component /5fc along the pi vector 
in the quark bulb is ~ 1 (similar to the case of upper quark bulb where the component oa: 
along the vector pi is ~ 1, see eq. (|36D). Therefore, ~ ^ and we see now that our 

/ 2 2 

integral over ap in the region w ^ ^ a' ~ — is indeed logarithmic. The lower limit 



>.ii^ ^^^^^^^ Y Cs ^ P 

of logarithmical integration is provided by the matrix element itself (since ~ 1 in the 
lower quark bulb) while the upper limit, at a^ ~ m? /C,s is enforced by the non-zero C and 
the result has the form 



iiTrij'^{pA_)ij^^{q^-pA_))) 



Fig. 6 



9"- 



d P'^ PI + p' [ 

4:11'^ P\p'\ 



-I^'ip'^ 



(57) 



where 



^In the region we are investigating, we can neglect the /3' dependence of the lower quark loop. 
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/^(j>1) = ^/^$^(«>i+p1) (58) 

is the impact factor for the lower quark bulb. It is easy to demonstrate that can be 
reduced to the double- integral form (|37| ) (with the trivial change pa Vb)- 

Let us compare now the matrix element ( |57D with the corresponding contribution to 
physical amplitude shown in Fig. ^ which has the form: 

~2 J 16^16^ 

{ap(3pS-p]_ + ieY{a'p(3'pS - p'l + ieY[{ap - a'p){(3p - (5'p)s - {p - p'Y + ie] ' 
where the upper quark bulb is the same as in Sec. 

refsec:Gen.structb. This integral is rather similar to the one for the matrix element of the 
operator, except that there is now a factor of the upper quark bulb, and there are integrals 
over Pp and p^. 

The previous arguments show that the logarithmic contribution comes from the region 
a;p^ap~^, ^~/3-C/3'-Cl, and now the upper limit of the logarithmic integral is set 
not by the regularized path-ordered exponential, but by the upper quark bulb, at ~ 1. 
Hence we have 

L.h.. of E,. (59) ^ >|! 1„ (^) /^;Mli±Pl,.(,J,B(,- ). (eo) 

2 

This agrees with the with estimate Eq. (|57|) , if we set = This corresponds to making 
the line in the path-ordered exponential have a finite rapidity relative to the photon. 

Thus, a more correct version of the factorization formula Eq. (|4l|) or Eq. (^4]) has the 

operators tj and f/t "regularized" at C ~ 

j d^x J d''z6{z.)e'P^-^T{j^{x + z)Uz)} 

= Ee?/ ^/i(Px)Tr{f/^-^(p)m^-^(-p)} 

i 

-|-term with higher order impact factors. (61) 



D. Next-to- leading order 

Let us outline the situation in the next-to-leading order in coupling constant. At the level 
of the 0{g^) calculation that we are examining, our results so far show that this formula 
captures all of the contributions from the original graphs (as s — > oo) except for those from 
Op ~ 1. As indicated earlier, in Sect. [Ill A] , these are included if we define the O^g"^) impact 



23 



factor suitably. The reason why we have indicated the higher order impact factor separately 
in Eq. ( pT] ) is that it is more than a trivial higher order correction to the lowest order impact 
factor J^^,, as we will now show. 

Following the strategy indicated by Eq. (|77|), we examine the difference 

J d'x J d'z6{z.)e'P^-^T{j^{x + z)Uz)} 

-Ee?/ ^/.1.(l'x)Tr{?7C-^(p)f/K-^(-p)} (62) 
in an external field. We are assuming a calculation to 0{g^) in Eq. (|62D . 




FIG. 8. Impact factor in next-to-leading order in Og. 

Typical diagrams are shown in Fig. ^ From the shock-wave picture of external field (see 
Appendix 0) it is clear that the general form of the answer for Eq. ( p^ ) is 

g^l d^x^d^y±d''z^J^{x^,y^,z±) Tr{rf/(xx + z±)t'U^ {z^)}[U{y^ + z^)U 

J d'^x±d^z±J^{x±) Tt{U{x± + z±)U\y±)}. (63) 

where [U]ab is the Wilson-hne gauge factor ( p3D in the adjoint (gluon) representation. The 
first term corresponds to the case when the shock wave hits two quarks and a gluon and 
the second to when it only hits two quarks. Without the subtraction term in Eq. (0), the 
diagrams in Fig. would diverge logarithmically at small ap. But after the subtraction the 
result will converge; the integral (to leading power) will be dominated by a ~ 1. 

Thus we obtain the coefficient functions (impact factors) Ji and J2. So, the operator 
expansion up to the next-to-leading term has the form 

J d'x J d'z6{z.)e'P^-^T{j^{x + z)Mz)})a 

= ^6^1 / d'z^I^,{x^)TT{U<{x^ + z^)U^<{z^)} 
+ g^J d^x±d^y±d''z^jf{x^,y±) Tr{t"f/^(a;x + ^±)t'f/^H^±)}[f>Hz/± + ^±)U 
+g'J d^x^d'z^4{x^)TT{U'^{x^ + z^)U^^{y^)} 

+0{g'), (64) 
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where the operators U in the 0{g'^) term must be also regularized at C = ^ in order to 
simulate a proper cutoff for the logarithms ~ gf^ln^ as well. In principle, this procedure 
may be repeated many times yielding the coefficient functions (impact factors) in any given 
order of perturbation theory just as for the usual Wilson expansion. 

It is worth noting that the above procedure of separating the Feynman integrals into the 
contributions coming from large and small components of the momentum pA can be repeated 
for the bottom part of the diagram with the result being the separation of loop integrals 
into contributions of large and small components along the pb- In the leading order in as 
the result will have the same form as Eq. (^Tj): 

= E / ^C(P^^ Tr{f/'(p)^>(g - p)}. (65) 

Here 

U'{x±) = [OOPB + X±, -OOPB + X±] , f/'t(x_L) = [-OOpB + X±, OOpB + X±] (66) 

are gauge factors ordered along a straight line approximately in the direction of motion of 
the lower quarks (pb) and the impact factor will be given by the same expression (|37D save 
the trivial change p\ p^. Therefore, the amplitude of scattering of virtual photons at 
high energy (|) can be represented as a product of two impact factors times the vacuum 
expectation value of four Wilson line operators representing the gluon ladders (see Fig. 

^rr'6^'\q^)A{pA,PB,q) = -^^(Z^f j ^^/^(Px, gx)/^(pl, 0) 

(0|Tr{f/(px)m(g^-p^)}Tr{{7'(pl)?7't(-pl)}|0). (67) 

The operators JJ are "normalized" in such a way that they are ordered along the slightly 
non-light-like line collinear to pa and U' along a line collinear to pb- 



IV. ONE-LOOP EVOLUTION OF EIKONAL-LINE OPERATORS 

In previous sections we demonstrated that at large energies the scattering amplitude of 
virtual photons can be reduced to the matrix element (between the "virtual photon" states, 
see definition (PD) of the two Wilson-line operators defined on a near-light-like line collinear 
to Pa- Now we must study the dependence of these matrix elements on energy which reveals 
itself through the dependence on the slope of the supporting line. So, we must find^ 



''In this Section we omit the trivial gauge end factors (Bq) which will be restored in next Section 
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FIG. 9. Factorization of gluon ladders from the high-energy ampUtude of virtual photon scat- 
tering. 

X—Tr{U{x^)U\y^)} = -2(—TT{U{x^)U\y^)} (68) 

at large A. The operators IJ and f/^ are defined on the lines coUinear to p'' = pi + C,p2 where 
C = ^ is a small parameter which determines the deviation of the supporting line from the 
light cone. This derivative can be expressed as 

C|Tr{f/(a;x)m(yx)} = 
igCJudu (Tr{[oo, u]^F^,{up'^ + Xj_)[u, -oo]^.f/1'(?/_L)} 

- TT{U{x±)igC Judu[-oo, u]yF^,{up^ + yi_)[u, oo\y]^ (69) 

where = [up^ + x±,vp'^ + x±]. So,the derivative of the two Wilson-line operator 

has reduced to a more complicated operator and therefore, in general, we can extract no 
information on the behavior with respect to A. But in the case of large A (small Q we can 
expand the complicated operator in r.h.s. of eq. (|69|) in inverse powers of A as it was done 
for the T-product of quark currents in the previous Section and we will show in this Section 
that the result will have a similar form 

igC, J udu 

(Tr{[oo,M]a;-F;,(Mp^ + x_l)[m, -oo]x.f/"^(y±)} - Ti{(j{x i_)[-oo, u]yF^,{up^ + ?/±)[m, oo\y]^ = 
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- 4^ Jdzs. 

(Tr{[/^(xx)f/tf(z^)}Tr{f/«(z^)f/t<(y^)} - N,TT{U<{x^)U^'^{y^)}^ ^^^^ 



+ Oig') + 0(C) (70) 

So, at large energies (A) the derivative of the Wilson-hne operator do reduce to another 
operator constructed from Wilson lines. Unfortunately, the number of Wilson lines is not 
conserved since the equation is non-linear. However, we shall see in Sect. 5 that in some 
important cases the eq. (|70D reduces to linear BFKL equation. In the rest of this Section 
we shall derive this equation which is one of the main results of this paper. 

In order to establish the operator equation (|70D let us compare matrix elements of the 
l.h.s and r.h.s. If we knew that the expansion goes in terms of Wilson lines beforehand, it 
would be enough to compare these matrix elements between two (or four) real gluons. But 
since we want to prove that the gluon operators in r.h.s of eq. (^0]) assemble in Wilson lines 
we must compare these matrix elements between an arbitrary number of real gluons - i.e., 
in external gluon field (see the discussion in Sect. 3). So, we must find the matrix element 
of the operator in l.h.s. of eq. (ffOD in the external gluon field at large A (small Q: 

igC!udu{Ti{[oo,uyj:F^,{up'^ + x±)[u, -oo]^.?71'(?/_l)} - 
TT{U{x^)ig( Judu[-oo,u]yF^,{up'^ + y±)[u, oo]y})A (71) 

In the lowest order in coupling constant we obtain zero since C ^ and the external field is 
independent of (. But already in the first order in the limit of the matrix element ([7T|) 
is nonvanishing due to the longutudinal divergencies. (As we shown in previous Section, 
some of the contributions to the matrix element of the operator TrU{x±)W{y±) contain In^ 
which means that the derivative is ~ ^ so the r.h.s. of eq. (|69D is actually nonvanishing 
at C ^ 0). Let us calculate the matrix element (^) in the one-loop approximation. It is 
convenient to use the light-hke gauge A^, = with the vector n directed along pb (although 
all the calculations can be repeated in the background-Feynman gauge with the same results, 
since we have checked that the contributions due to gauge terms ~ in the propagator in 
the axial gauge cancel). 



A. Calculation of the diagram in Fig. 10a 

In the first order in there are two one-loop diagrams for the matrix element of operator 
(|69|) in external field (see Fig. p^O]) . 



We shall start with the diagram shown in Fig. |l^a. The calculation is quite similar 
to the calculation of impact factor considered in previous Section. For the sake of future 
applications we shall calculate the derivative C-§^U{x±) tJ^{y±) where tJ{x±) tJ^{y±) = 

{U{x±)yj{U^{y±)}^ is a product of Wilson-line operators with the non-convoluted color 
indices. First, using the expression for the axial-gauge gluon propagator in the external field 
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A B C 

FIG. 10. One-loop diagrams for the evolution of the two- Wilson-line operator. 



from Appendix C we obtain? 

igC Judu{[oo,u]^F^,{up'^ + x±)[u, -oo].^ (g) f/"^(i/±))A 
= gC J udu J du 

Q{u' — u) [oo, ujxf^lu, u]xF^:,{up^ + x±)[u, — oo] 
® / dv[—oo,v]yt^[v,oo]y 



(up< + ipj - V.^)0^\p^, - ^V.) 

VV ^ P-P2 P- P2 



vp^ + yi. 



ah 



+ gC du[oo, ulxt"" [u, —oo]x ® / dv[—oo,v]yt [v,oo] 



VV « p.p^ ' P-P2 



vp^ + yi. 



ab 



(72) 



where the operator O has the form: 



O =4—F^—F,- -(D'^F + 



^D-F^^ - J^V^D'^F^sJ^)^^ 

(73) 



P-P2 P- P2 2p ■ P2 2p ■ P2 V'^ 



(With our accuracy, multiplication by coincides with multiplication by pi). We may drop 
the terms proportional to V, in the parenthesis since they lead to the terms proportional to 
the integrals of total derivatives, namely 



/ du[oo, uj^lu, -oo]p^(L'^<I>(mp^, ...))ab = 

J du£{[00,u]t-[u,-00]mup^,...))ab} = 



(74) 



^It can be demonstrated that further terms in expansion in powers of gluon propagator (C5) beyond those 
given in eq. (|C^ ) do not contribute in the hmit C ^ 
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and similarly for the total derivative with respect to v. Therefore we may rewrite eq.(|72D as: 
= gC J udu J du' 

+ Q{u' - u)[oo, ujxf^lu , u]xF^,{up'' + x±)[u, -oo] 



(g) / dv[—oo,v]yt [v,oo]i 



vp'^ + yi. 



+ gC J du[oo,u\xt"'[u, —oo]x ^ J dv[—oo,v]yt^[v, 
i^up'' + x± p^O, 



oo 



vp^ + yi. 



(75) 



Now let us consider the limit C ^ 0. The Wilson lines made from external fields are 
regular in this limit and the only singularity which can compensate C, in the numerator of 
l.h.s of eq. ( |75D is ^ coming from the differentiation the gluon propagator in the external field 
which contains terms ~ In C,. Therefore only the last term in eq. (|75| ) gives the nonvanishing 
result . Adding the similar contribution from the second term in r.h.s. of eq . (|69|) we have: 

C^^{ij{x^)U\y^))A 



(76) 



g I du[oo,u]j:"'[u,—oo]x® j dv\—(X),v\yt [v,oo\ 



{upA + X_L Up^O,, - vO,,p^ 



vpA + y± 



ab 



Let us at first neglect the gauge factors [oo, ul^lu, — oo] and [— oo, , oo] - in other words, 
let us consider the trivial zero-order term of expansion of these gauge factors in external field. 
We have then 



^g'l 



b(2\2 



[ak5\(3k + C«fc)5(A + CaO + a;5'(A + C^dm + Cafc)]e-^(''")^+*(''^)^ ((fc|a.|/))^^ (77) 

Note that as in the case of quark propagator we need the Green functions integrated over 
the a component of external field. The calculation of fast-moving gluon propagator in the 
external field mainly repeats the derivation of the formula (|3^) for the quark propagator and 
we will only sketch it here. In the lowest order in expansion of the operator ^ in powers of 
in eq. ( [75D one obtains: 



(7J 



{k 0„k-p] 

= -TTT- \-D''F^,{p) + 2 / (pO ]_ ■ F^.iP - P) 

ttk + le [ J IGtt^ —jj'p + leak 



{k — pY + ie 
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As can be seen from eq. (^) our j3p are ~ C«p so we can neglect them in the arguments of 
external field. Then the expression in braces is proportional to one of the operators: 

[DF]{x^) + 2i[FF]{x^) 

= Jdu[oo,u\xD°'Fa,{upi + x±)[u,-oo]x + 

2i j du jdvQ{u — v)[(yo,u]xF^{upi + xj^)[u,v]xF^,{vpi + xj^)[v,—(X)]x (79) 
at afc > or 

[DF^]{xi_) -2i[FF^]{xi_) 

= Jdu[-oo,u]xD"'Fa,{upi+x±)[u,oo]x + 

2i j du j dvQ{v — u)[—(X),u]xF^{upi + XA_)[u^v]xF^,{vpi + x:[_)[v,oo\x (80) 

at Ofc < in the lowest order in external field. It can be demonstrated that the subsequent 
terms of the expansion of the operators in eq.([75|) in powers of external field "dress" 

lowest-order expressions for [DF] and [FF] by proper gauge factors according to eq. ( [75|J5(]| ) 
so we have 



s -p j~\ rp s \ I 



k 



p 



ab 



2iTi5{ap) 



Qiak){i[DF]{p^) - 2[FF]{p^)} + ei-ak){^[DF^{p^) + 2[FF^{p^)} 



1 



(fc-p)2+ze 



It can be simplified even more if one notes that the operators in braces are in fact the total 
derivatives of U and W with respect to translations in the perpendicular directions: 

9lUix^) = 9&-Uix^) = -t[DF]ix^)+2[FF]ix^) 
dlUix^) ^ sA^W{x^) = t[DF]{x^) + 2[FF](xx) 



(82) 

(note that dj_U = —d'^U). So, the eq. ( pT|) reduces to the expression for the fast-moving 
gluon propagator in the external field in the form: 

(, a. k-p))= -2"^K)9(«0ait/(Px)-9(-..)aic/t(pj 

^-^ ak [k'^ + ie)[{k — p)'^ + ie\ 

This formula is valid in the region (|2l|) provided the f3 component of the overall momentum 
transfer to external field is small (in our case Pp (, see eq. (^)). Substituting now eq. (|83|) 
in eq. (^) we have: 

„2^a ^fbr (fik^ (fil^ ^-i(k,x) , +i(l,y) f doj^ 
J 47T-2 Att^ j Ott 



^ k A 



{(als+kl)^ 



Ca^s+kAQic.)id'U{k^ - h)) - e{-a){d'U\k^ - /x)) (^iflf^ 

^ fc _L fc _L'' 



1 

-1,2 c-i_i.2 



r dk^ dp± 1 
47r J 47r2 4it^ k'^ (k-p)'^. 



[{d^U{p^))al> + {d^U\p^))ab] 



,-i{k,x-y)±~i{p,y)± 
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where the integral over converges at ~ ~ 1. 

Now let us turn our attention to the omitted gauge factors [oo, u\t"-[u, —oo] and 
[oo,f]t*[f, — oo] in our starting expression ([751). We demonstrate in Appendix |B| that they 
should be substituted by )f:"[oo, — oo] ® t^[—oo, oo] or [oo, — oo]t" ® [—00, oo\t^ depending on 
the sign of a^- (In the coordinate space it means that the transition through the shock-wave 
" wall" can be before or after emission of quantum gluon depending on the sign of x^, and , 
see Appendix After that our final result for the contribution of the diagram in Fig. p!o| a 
reads 

C^(f/^(x^)m^(l/x))A = 



3_ 

in 



(( 



x± 



^5) 



Using the identity 

it can be rewritten in the form: 



^6) 



2: 
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{Wiz^)Uix^)}^{Uiz^)Wiy^)}l + {Uix^)Wiz^)}l{Wiy^)Uiz^)} 
S^{Uix^)Wiy^)}l + Si{W{y^)Uix^)}^ 



[x-zy-z)^ 
{x-zy-' (y-z)'^ 



i7) 



where we have displayed the color indices explicitly. 



B. Calculation of the diagram in Fig. 10b 



The contribution of the diagram in Fig. pll b is calculated in a similar way. One starts 
with the expression (cf. eq.(|73|)): 

= gC J udu J du 

{Q{v - v')Q{v' - u)[oo,v\xt"-[v,v'\j:t\v' ,u]xF^,{up^ + -oo\x 
+ Q){v - v')Q{u - v)[oo, u]^F^,{up'^ + xj_)[u,v]J''[v,v%t''[v' , -oo]^ 
+ Q{v — u)Q{u — v') [00, vjxf^lv, u\xF^:,{up'' + x±)[u, v']xt^[v', —00] 

®U\y±) 
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vv ^ r) ■ r)r, ' n ■ rir. 



P2ri 



P-P2 



P-P2 



vp^ + X± 



ah 



+ 



ig(, J udu[oo, u]xF^,{up'^ + x±)[u, — oo]^; (g) 
dvdv'Q{v - v')[oo, vjyf'lv, v']yt^[v', -oo] 



P2v 



(m/ + {pi - V.-^)0^'>{j)^ 

VV ^ P-P2 P-P2 



V.) 



Vp^ + X_L 



ah 



+ 



ig'^ J dudv {Q{u — f )[oo, u\J,"'[u, v\xt^[v, — oo]^; 



+ x^ pM - - -^v., 

\ ^ P-P2 P-P2 



vp^ + X±_ 



ah 



(88) 



First, let us demonstrate that the contribution of the terms ~ V, in the parentheses in r.h.s. 
of eq.(|88D vanishes (cf. eq.(|73[)). Indeed, using eg. ([7^) and integrating by parts it is easy to 
reduce these terms to the sum of the contributions of the type 



g'^fabcl-oo du[00, ulJ^'lu, -OO]^ 



+2iFi.^ V^g^x+2iFi.^ 



P-P2 P-P2 1^ P-P2-' V^gpri + 2lFprj } p-p2 

This expression corresponds to diagram shown in Fig. |TT[ 



vp'^ + x± 



ah 



(89) 




FIG. 11. Typical loop integral for the contribution of the gauge terms ~ n \\ P2- 
Let us consider the integration over (3k in the loop integral corresponding to the mo- 



mentum k. As we shall see below, similarly to the case of diagram shown in Fig. p!0| a the 
contribution which survives in the limit C ~^ has characteristic ~ 1 while all other a' s 
corresponding to external field are ~ ^ ~ C- This means that in the contour integral over 
(3k all poles coming from the denominators 
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{{ak + a){(3k + f3)s-{k+p)l + ie}-' 



(90) 



lie to the one side of X axis so that the resuhing integral is zero. This happened since we 
have cancelled the eikonal denominator in the original diagram in Fig. p!0|b. We shall see in a 
minute that this eikonal pole in f3k may lie to the opposite side of the X axis thereby leading 
to the non-zero contribution for the terms not proportional to V,. Actually, the cancellation 
of the terms proportional to longitudinal part of the gluon propagator (in the external field) 
is a consequence of the gauge invariance of the operator U. 

So, we have reduced the contribution of the diagram in Fig. pil| b to 



= gC J udu J du 

{Q{v - v')<S){v' - u)[oo,v\J''[v,v']xt\v' ,u\xF^,{up'^ + xa_)[u, -oo]a 
+ Q{v - v')Q{u - v)[oo, u]:,:F^,{up'^ + X _i_)[u, v]j"-[v , v%t''[v' , -oo]x 
+ {Q{v - u)Q{u - v')[oo, v] J^'lv, u]^F^,{up^ + x±)[u,v%t^[v', -oo]^ 

^up'' + x± 0„ vp^ + 

+ 

igC, j udu[oc, u\xF^, + {up''- + x_l) [m, — oo]^; ® 
dvdv'Q{v - v')[oo,v]yt"-[v,v']yt''[v' , -oo]y 



x± 



Vp^ + X_L 



ah 



+ 



ig"^ / dudv {Q{u — v)[oo,u\xt'^[u, v\xt^[v, —oo\a 



(^Up'^ + X_L p^O, 



Vp^ + X_L 



(91) 



As in the case of the diagram in Fig |10pL, a nonvanishing result comes only from the last 
term where we differentiate propagator in the external field thus obtaining an extra factor 
= 1 in comparison to ^In^ — for the first three terms. We have then 



igC J udu{[oo,u]a:F^,{up^ + x±)[u, -oo]x)a = 

/oo ru 
du \ dv\oo, u\xt°'\u, v\x^\v, —'x^xVu.pA + x^ up^O,, — vO,,p^: 
-oo J —oo ^ 



VPA + X_L 



ab 



(92) 



where we have omitted ®U'^{yx) for brevity (it will be restored in the final answer). 

Again, let us neglect at first the gauge factors [oo,m], and [f, —oc\. We have then: 
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ab 



(93) 



Now we can use our result for gluon propagator 



dak dPk dk± dp± 



-i(p,a;)j 



and obtain 
1 



{Pk +^ Oik - if-Y akPkS -k\ + ie 
[e{a){d^U{p^)) - Q{-a){d^U\p^))U- ^ 



An 
47r 



dkx dpj 



471-2 4^2 

1 
p 



-dlU{pj 



1 



{k-pf_ 



1 

pi 



ab 



(94) 



In the Appendix ^ it is demonstrated that the gauge factors of the type [m, v\ which we 
omitted lead to the substitution [oo, u\t"-[u, v\t^[v, —oo] — >• f^U {x±)t^ (in the coordinate space 
it means that the interaction with the shock-wave occurs between emission and absorbtion 
of the quantum gluon). So we finally obtain: 



9^ 
-—t 
An 



dz I 



p^ p^ 

U{z±) Tt{U{x^)U\z^)} - N,U{x±)] ® U\y±) 



(95) 



{x - z)l 

where we have restored the omitted factor W{y±). The contribution of the remaining dia- 
gram in Fig. pUj c differs from eq.(P3|) only in the substitution U ^ 



-.2 
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dz^U{x^) ® U^{z±) Tr{f/(zx)f/t(y^)} - N,U^{yj 



(96) 



Now we are in a position to write down the final answer for the one-loop evolution of the 
operator U{x±)U^ {y±). Combining the expressions (p2D, (pS]), and (p^) we obtain: 

= Jdz^{[{UHz^)U{x^)}'^{U{z^)UHyM + {U{x^)UHz^)}l{UHy^)U{z^)}'; 

,fcl {x-z,y-z)^ 



-6^{Uix^)U\y^)}l - 6l{U\y^)Uix^)y^ 



{Uiz^}}) Tt{U{x^)U\z^)} - N4U{x^)y^ U\ 



[y±n 



{U{x^)}) \u\z^)fTT{U{z^)U\y^)} - N,{U\y^)Yi 



(x - z)l 

1 

(y - z)[ 



(97) 
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This is the one-loop result for the operator U'^'^{x±) ® f^^'^^(y±) in the low-a external field 
A corresponding to the bottom part of the diagram in Fig. 2a. As we discussed in previous 
Section, the operator form of this one-loop evolution is: 

C-^{uHx±)}]{u^Hy±)}1 



1 



(x - z)l 

1 



(98) 



where the operators f/^ and f/"!"^ are integrated along the line coUinear to in order to 
impose cutoff a < in the matrix elements of these operators. 

V. EVOLUTION OF EIKONAL-LINE OPERATORS IN LEADING LOG 

APPROXIMATION 



A. Linear evolution at large 

Let us outline how to obtain the energy dependence of the amplitude using the expansion 
in eikonal-line operators. As we have discussed in previous Sections after formal expansion 
at large energy we obtain in the leading logarithmic approximation the operators U and U'^ 
"normalized" at the slope C = ~ times the impact factor: 



JdxJdzS{z.)e'P^^T{j,{x + z)Uz)} = e]J dx^dzj^,{x^) Tr{f/^(xx + z^)U^^{z^)} + 0{g' 

(99) 



flavors 



where /"^(xj,) ig given by eq. (40) and the dots stand for the next-to- leading term given in 
eq. (|6^) . (Hereafter we will wipe the label {") from the notation of the operators). Matrix 
element of this operator -C U^{x^)U^^{ii^) ^ (see eq. (|5lD for the definition) describes the 
gluon-photon scattering at large energies ~ s. The behavior of this matrix element with 
energy is determined by the dependence on the "normalization point" C,. From the one-loop 
results for the evolution of the operators U and f/^ it is easy to obtain evolution equation: 
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C-^TT{U^ix^)[x^,y^].U^<{y^)[y^,x^] + } 

,2 



9 



iV,Tr{f/^(x^)K,yJ_f/tC(y^)[y^,a;J+})- (100) 



where we have displayed the end gauge factors ( pSf ) exphcitly. We see that as a result of 
the evolution the evolution of the two-line operator Tt{UW} is the same operator (times 
the kernel) plus the four-line operator Trjf/f/"''} Tr{?7f/^}. The result of the evolution of the 
four-line operator will be the same operator times some kernel plus the six-line operator of 
the type Tt{UW}Tt{UW}Tt{UW} + Tt{UWUW}Tt{UW} and so on. Therefore it is 
instructive to consider at first the large- A'^c case where, as we demonstrate below, the number 
of operators U is always the same during the evolution. 

Let us consider the evolution equation ( |69D in the large- iVc limit. It is convenient to 
rewrite it in terms of the operators 

V{x^,y^) ^ ^^{x -y)l'iTT{Uix^)[x^,y^],U\y^)[y^,x^]+} - N,) (101) 

so it reduces to: 
d 

C-g^V{x±,y±) 

O-N^ !^^^ [pElllA + ZlfilJ^ _ n^x,i/x)(^x-yx); ^ ,102) 



167r3 ] y±y {x± - z±y {x± - z±y{z± - y±)'^ 

It is easy to see that matrix elements of the operator ((V^^)) are ~ in comparison to the 
matrix elements of the operator {{V)) ! 

So, in the leading order in A^^ the evolution of the forward matrix element of the two-line 
operator V{x±,y±) is governed by linear BFKL equation ^ 

4r,.)),^-£,./.4M^.^-i^) ,.03) 



'The only exception is disconnected contributions of the type {{V)) (g) (0|F|0) but they are 0{g^) corrections 
to the matrix elements {{V)) which we shall not consider in the leading logarithmic approximation. In the 
next order in they must be taken into account together with 0{g^) corrections to impact factor ( |40| ) 

^ The connection between Wilson-line operators and the BFKL equation was first discussed in Ref.El . 
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where {{V{x±))) = {{V{x±,0))), see eq. P^). The eigenfunctions of this equation are powers 
(a;j_)~2+*^ and the eigenvalues are —^N^xi^), where x{^) = ~ Re^/'(| + iz/) — C. Therefore, 
the evolution of the operator V takes the form: 

{{V<^{x^))) = I^ixl)-'^^- (S)'"''^'^'^ Jdz^{zlr'^-{{V^^{z))) (104) 

We may proceed with this evolution as long as the upper limit of our logarithmic integrals 
over a which much larger than the lower limit — which is determined by the lower 

quark bulb, see the discussion in Sect. 3. (In other words, if we look on the derivation of the 
evolution equation given in previuos Section we can see that it holds true as long as we can 

at a certain point Co such as 



neglect ~ — in comparison to ~ in eq. (p4D). It is convenient to stop evolution 



Co = A> < 1, g^lna < 1 (105) 

Then the relative energy between the Wilson-line operator V^^ and lower virtual photon will 
be So = m^cT^ which is big enough to apply our usual high-energy approximations (such as 
pure gluon exchange and substitution g^i, j^P2nPiu) but small in a sence that one does 
not need take into account the difference between gf^ln^ and o^ln^-j. Then finally the 
evolution ( |103|) takes the form: 

{{V<=^ix^))) = /d.^(zi)-^-((y^"(.J)) (106) 

i/ ^ Ji lit J 

Now let us rewrite this evolution in terms of original operators UU'^ in the momentum 
representation. One has then: 

((Tr{[/^=^(p^)[/tC=^(-p^)})) 

where we omit for brevity the end factors (|35|). Since we neglect the logarithmic corrections 
~ ^f^lncr matrix element of our operator U''°W'^° coincide with impact factor up to 0{g'^) 
corrections: 

((Tr{f/^°(p^)f/t^o(-p^)})) 

jiV^^ _ 1 2 r da ^^{ttppi - CoapP2 + P±) 



,2" 



2 ^^'J ns iCoal+piy 



= 9'^^T.el^l'{p^) (108) 
Combining eqs. (0)J^1O^), and ( |108|) we reproduce the usual leading logarithmic result 



for virtual 77 scatterinj 
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A{Pa,Pb) 



At s — s> the amplitude ( |109|) is determined by the rightmost singularity in the v plane 
located at z/ = (in terms of complex momenta plane it corresponds to the position of the 
"bare pomeron" at j = 1 + ^Nchi2) and the answer is 

(110) 

where C(3) — 1.202. In the case of small-x deep inelastic scattering the evolution of the 
matrix element (A^|V|A^)is the same as eq. ( |109| ) with the only difference that the lower 
impact factor should be substituted by the nucleon impact factor determined by the 
matrix element of the operator UU'^ between the nucleon states f 



(iV,p^|Tr{f/^°(a;^)t/t^"(0)}|iV,p^ + /3p2) = 2i,5{(3) / -^e^(^'-)^_/^(p^) (111) 



dp±^i(px)^ 1 

where 27rS{f3) reflects the fact that matrix element of the operator containes unrestricted 
integration along p'^°, (cf. eq. (^21) ). The nucleon impact factor I^{p±) defined in (|111|) is a 
phenomenological low-energy characteristic of the nucleon. In the BFKL evolution it plays 
the role similar to that of nucleon structure function at low normalization point for GLAP 
evulution. In principle, it can be estimated using QCD sum rules or phenomenological 
models of nucleon. 

Let us discuss how the nucleon impact factor ( [111| ) is related to the gluon structure func- 
tion of the nucleon which is defined as the matrix element of the gluon light-cone operator: 

V{uj,^i) = -— fdue-'^^^^{N\TT{F^{up,)[up,,0]F^.[0,up,]}\N)'' (112) 

SUl J 



'^This is called "hard pomeron" contribution to the structure functions of deep inelastic scattering since all 
the transverse momenta in our calculations are large or ^ which is the same in leading logarithmic 

approximation). However, due to the so-called diffusion in transverse momenta the characteristic size of the 
in the middle of gluon ladder is eV^^^^ (see e.g. Ref.El for discussion) so at very small x the region 
p± ~ ^QCD may become important. It corresponds to the contribution of the so-called "soft" or "old" 
pomeron which is constructed from non-pcrturbative giuons in our language and must be added to the 



hard-pomeron result given by eq. (110) 
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where /i is the normahzation point for the hght-cone operator. (The unrenormahzed operator 
F{upi)F(0) is UV divergent so we regularize it by counterterms just as for the local operator, 
see e.g. EJ ).The physical meaning of fi is the resolution in the transverse size of the gluon: 
T>g{uj,fi) is the probability to find inside a nucleon the gluon carrying the fraction a; of the 
nucleon momentum with the transverse size fi~^. Formally, 

ldk^Qifi^~kl)Vg{LU,k^) (113) 

where Vg{uj,k±) is the gluon distribution over transverse momentum k± and fraction of 
logitudinal momentum upi. 



ujVg{u},x^) = -- J due-'i^''{N\ Tt{F^{upi + x±) 

[upi + x^, , 0]F5.(0)[0, -oopi + x^]}\N) (114) 

It is easy to relate impact factor to the gluon distribution - actually, it is the same quantity 
with different regularization of longitudinal integrations. Indeed, 

2 J, JN\ TT{d,U<oip^)d,W<^i-p^)}\N) 

jj ^^^^ = 2 2^) 

4 f 

= — g"^ / du{N\ TiU-oop'^" + m/" + x±]F^{up'^" + x±) 
s J 

+ x±, -oo/o + x±][-oop'^", 0]F^,{0)[up'^'\ -oo/°]}|A^) (115) 



Now we see that the r.h.s.'s of eqs. ( [114| ) and ( |115|) coincide up to a different cutoff in the 



longitudinal intergation in matrix elements: in the case of gluon distribtuion the integrals 
over the a component are restricted from above by ^ whereas for the matrix element ( p.l5| ) 



the cutoff is a/ ?^ = — so they coincide at a = u . Therefore 



'For example, in the case of diagram in Fig. ^the contribution to impact factor (115) is (cf.eq. (55)) 

^dopda^dp^dp^^ g^N,{N^ ^ l)plT:,ip, -p')TUp, -p')<^gip') 

2tt 2tt 2tt 4tt'^ {(a'^s +p\- ie)^{a'pf3'pS - p'l +ie)^[{a - a'){l3 - P')s - {p - p')^ + ie] 



(116) 



whereas the contribution to the gluon distribution defined by r.h.s. of cq. (114) has the form: 



i f dap da'p dPp dp'^ 



A J 2n 2iT 27r Att^ {aojs — pj_ — ief 



, _ g^A^c(iv.^ - i)pir-.(p, -y)rg.(p, -pO$f,(y) 

^"^^ '^'\a'p(5'pS-p'l+ieY[{a-a'){(3-(3')s-{p-p'f + it\ ^ ' 



so we see that the only difference is in the cutoff for the logarithmical integration over a 
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^/^(p^)=aP,(a,p^) + 0(/) (118) 



pi 



where the impact factor is determined by Wilson hnes U and f/^ parallel to pi + -^crp2 1 1 
(yp2 + —Pi and airP' plays the role of the relative energy between Wilson lines and nucleon. 



B. General case 



Unlike the linear evolution at large Nc, the general picture without large- iVc simplification 
is very complicated: not only the number of operators U and f/^ increase after each evolution 
but they form more and more complicated structures like displayed (and not displayed) in 
eq. ( |122|) below. In the leading log approximation the evolution of the 2n-line operators 
such as Tr{t/t/^} Tr{?7?7^}... Tr{[/[/''^} come from either self-interaction diagrams or from 
the pair-interactions ones (see Fig. |12]) 





FIG. 12. Typical diagrams for the one-loop evolution of the n-line operator. 



which we have already calculated in Sect. ^\ — see eqs. (|85| , |95| , |96|) . Therefore the one- 



loop evolution equations for these operators can be constructed using these rules which we 
shall list here for completeness in the explicit form (in this Section we will use the notation 
Ux = U{x±) etc): 

d 



' dzA{um';{UMl}l + {u.uljiiulu.}'; - 6';{UxUl}i - sHulu.}';: 



[X 



z,y 



(X — z) 
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dz^imiiUyUiu.}^ + {u.uiUyjim'; - miiUy}^ - {Uyjim'r, 



dz^mVAuluM)] + {ulu^ulYM}] - m\{ul]] - {ul]\{ul] 



{X- z,y- z)± 
ix-z)l{y-z)l 

{x-z,y-z)± 
ix-z)l{y-z)l 
(119) 



for the pair-interaction diagrams in Fig. 12a and 



^ / dz^ PI Tt{U^WJ - NMl] j^^y^ 



(x-z)l 



(120) 



for the self-interaction diagrams of Fig. |T2|b type. Using eqs. (|119| , |120D it is easy to write 
down evolution equations for arbitrary n-line operator. For example, the evolution equation 
for the four-line operator appearing in the r.h.s. of eq. (|100|) has the form: 

Tr{f/,.[x, [z, x] + } Tr{f/,[z, z] + } 



|t/tx{[Tr{t/,.[x,t]_f//[t,x] + }Tr{f/t[t,^]_f/t[^,t] + }- Ar,Tr{?7,[x,^]_[/][^,x] + } 
{x - z)l 



IGvr^ 

i:i{U-Az,y]-Ul[y,z]+} 



■^-'--^(x-t)i(z-t)i 
+ Ti{a,\x,z\-Ul[z,x\ + } 

[Tr{f/4z, t]_f//[t, ^] + } Tr{f/4t, y]_f/t[y, + } - iV, Tr{f/,[z, ^] + } 

+ [Tr{[/,[x,z]_f/t[^,t]+f/4t,y]_t/t[t/,^]+f/,[^,y]_f/t[t,x] + }+ 
TT{U,[x,t].Ul[t,z]+U,[z,y].Ul[y,t]+Ut%z].Ul^ 

{x-t,y- t)_L 1 {x-t,z-t)± {z-t,y- 1)± 

- t)l{y - t)l {z - t)l + (x - t)liz - t)l ^ (z- t)liy - t)l 



{y-t)l{z-t)l 



'1211 



where we have displayed the end gauge factors (0) explicitly. Note that each of the separate 
contributions (|119|) and ( p.20| ) corresponding to the diagrams in Fig. [l^ a and |T2| b diverges 
at large z while the integral (over t) in the total answer (|7^) is convergent. This is a case 
of the usual cancellation of the IR divergent contributions between the emission of the real 
(Fig. 0a) and virtual (Fig. |l2|b) gluons from the colorless object (corresponding to the l.h.s. 
of eq (0)). Another example of this cancellation is eq. (|100|) , see the discussion above. 
So, the result of the evolution of the operator in the r.h.s. of eq.(^) looks like: 
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Tr{[/^.[x,i/]_[/t^[i/,x] + } 



oo 



^0 2vr^ Co J 

[A^{x,z\z',...z\y)TT{U^>[x,l]M^'[l,xU} 
Tr{Ut[l,2]M^°[2, 1]+} . . .Tr{[/^«[n,2/]_[/f'[y,n]+} 

Tr{f/i<' [x, l]_f/|f° [1, 2]+f/|° [2, 3]_f/f ' [3, l]+f/f« [1, 2]_f/t^« [2, 
Tr{[/3^°[3,4]_f/i^«[4,3]+}...Tr{f/^«Ky]_f/t^o[^,n] + } + ... 

+ iV,"a(a;, z\ z\ ...z\ y- ) Tr{[/i° [x, y]_f/f ' [y, (122) 
where U^^'^ = U^'^\z^), = and 

An{x, z^, z"^, ...z^, y), Bn{x, z^, z^, ...z", y), . . . C„(x, z'^, z"^, ...z"-, y) are the meromorphic func- 
tions that can be obtained by using the eqs.( p.l9| , p.2CI| ) n times which give us a sort of Feynman 

rules for calculation of these coefficient functions. If we now evolve our operators from C ~ 
to Co given by eq. ( |105| ) we shall obtain a series (|122|) of matrix elements of the operators 



(f/)"(?7^)" (see eq. (|75D ) normalized at Co- These matrix elements correspond to small energy 
~ m? and they can be calculated either perturbatively (in the case the "virtual photon" 
matrix element ) or using some model calculations such as QCD sum rules in the case of 
nucleon matrix element corresponding to small-x 7*p deep inelastic scattering . It should be 
mentioned that in the case of virtual photon scattering considered above we can alculate the 
matrix elements of operators UU'' ...UW perturbatively and therefore in the leading order in 
as we can replace all t/'s (and W's) except for two by 1. (Recall that U = 1+ig J A^dx^ + ... 
so each extra U — 1 brings at least 0{g)). So, we return to the BFKL picture describing the 
evolution of the two operators UU'^ similarly to the large N^. case. The non-linear equation 
( |1(JU ) enters the game in the situation like small-x deep inelastic scattering from a nucleon 



when the matrix elements of the operators f/f/^...f/f/^ are non-perturbative so there is no 
reason (apart from large N,^ to expect that extra U and U'^ will lead to extra smallness. In 
this case, at the low "normalization point" Co one must take into account the whole series 



of the operators in the r.h.s. of eq. (|122|) which means that we need all the coefficients 
a„, hn...Cn which must be obtained using the non-linear evolution equations (I100| , |104| ) etc. 

The situation may be simplified using Mueller's dipole pictures . Tech- 
nically, it arises when in each order in asln(^) we keep only the term 

Tr{[/f«f/|^°}Tr{t/i^°f/|^"}...Tr{[/fo[/tCo} - subtractions "In r.h.s. of eq. (|75|) - for exam- 
ple, in eq.(|7^ we keep the two first terms and disregard the third one. In other words, we 



take into account only those diagrams in Fig. ^ which connect the Wilson lines belonging to 
the same Tr{t/fc[/2+i}- (This corresponds to the virtual photon wavefunction in the large- A'c 
approximation). Then the diagrams of the corresponding effective theory are obtained by 



"By "subtractions" we mean this operator with some of the Tr{t/fc[/2+i} substituted by 



42 



multiple iteration of eg. (|69D and give a picture where each "dipole" Trjf/fef/^+i} can create 
two dipoles according to eg. (|69|) . The motivation of this approximation is given in Ref.0 . 



VI. CONCLUSION 

Let us summarize our results for the operator expansion for high-energy scattering. It 
has the form: 



2 

a„(x, z\ z\ ...z", y- In ^) Tr{f/'^,[a:, l]-f/^[l, + } 



dx J dz6{z,)T{jf,{x + z)j,y{z)} 

oo „ 
n=0 

Tr{[/^i[l, 2]„f/l^[2, 1] + } . . . Tr{f/^„[n, y].Ul<[y, n]+} 



2 

+ bnix, z\ z\ ...z", y- In ^) Trlf/'^^.fx, l]-f/^[l, 2]+f/^2[2, 3]_ 



t/l^[3,l]+f/^i[l,2]„f/l^[2,a;] + } 
Tr{t/f3[3,4]_?7i^[4,3]+}...Tr{[/^„Ki/]_[/tqy,n]+} 



2 

m 



N:cn{x,z\z\...z\yM—)'^^{U'^.[xM-Ul%,x]+] 



(123) 



where the notations are the same as in eg. (|122|) . The coefficient functions a„, 6„, ..c„ absorb 
all the information about the high-energy (A oo) behavior of the amplitude while the 
matrix elements of the Wilson-line operators, however complicated, are low-energy hadron 
characteristics. In terms of functional integral representation for the amplitude (0) we 
make a decomposition of all the fields into large rapidity fields (with light-cone fractions 

a > \J^) and small-rapidity ones (with a < -y/^)- The integration over large Sudakov 

variables {a > -y/^) gives us the coefficient functions (a„...c„)(x_Lfc) while the integrals over 

small a < form matrix elements of Wilson-line operators. Coefficient functions contain 
logarithms of energy In while matrix elements - only In (. The dependence on ( cancels 
in the final result and In ^ emerges just as in the case of usual Wilson expansion where the 
dependence of the coefficient functions and matrix elements on the normalization point is 
cancelled in a similar way: In % + In ^ = In % . 

M P P 

In order to find a dependence of the amplitude on energy using this operator product 
expansion we must proceed as follows: first, we integrate over light-cone fractions a ~ 1 - it 
gives us the operator Tt{UW} "normalized at the slope C = Second, using the evolutions 
eguation we reduce the two Wilson-line operators collinear to pa to the sum of the many- 
Wilson-line operators (almost) collinear to pb times the coefficient functions containing In 
In the leading logarithmic approximation, the evolution egs. ( |119| , |120D are enough; beyond 
that, one must consider higher-order corrections to these evolution eauations. Finally, we 
must compute the matrix elements of many- Wilson-line operators sandwiched between our 
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target states. In perturbation theory (e.g. for the deep inelastic scattering from the virtual 
photon) or at large A^^^ the evolution of the two- Wilson-line operator is enough - others have 
the matrix elements smaller by g'^ (or N(^- and we have the linear BFKL evolution described 
by eqs. ( |103| , |104| ). If the target states are non-perturbative (e.g. nucleons for deep inelastic 
scattering at small x) we must take into account the whole non-linear evolution ( p.22| ) even 
in leading logarithmic approximation. 

There is, however, one important difference between operator product expansion for deep 
inelastic scattering and our expansion for high-energy scattering. In the case of Wilson's 
expansion the coefficient functions were purely perturbative (up to possible contributions 
from small-size vacuum fluctuations c3 ) whereas all the non-perturbative dynamics was 
hidden in the matrix elements. This is not the case for our operator product expansion - 
both coefficient functions and matrix elements can have perturbative and non-perturbative 
terms. For Wilson's operator product expansion this perturbative vs. non-perturbative 
separation was due to the fact that it corresponds to the separation of the integrals over the 
transverse momenta: p\ > /i^ form coefficient functions and < fi"^ matrix elements ( and 
the characteristic scale of the coupling constant depends on the scale of transverse momenta). 
For the same reason, separation of integrals over longitudinal variables has nothing to do with 
scale of as - it is determined by scale of p± which can be either lagrge or small independent 
on longitudinal momentum. So, since both matrix elements and coefficient functions can 
have the contributions from small and large momenta - both of them do have perturbative 
and non-perturbative parts. We have of course calculated only perturbative contribution 
toi the coefficient functions which comes from the region of large p± - the non-perturbative 
contribution comes from p_i_ ~ Aqcd and it corresponds to the soft-pomeron contribution 
to the coefficient functions. So, in order to separate perturbative physics from the non- 
perturbative one for the high-energy scattering we must do the additional job and split the 
integrals over the transverse momenta in hard and soft parts both in coefficient functions 
and matrix elements. The study is in progress. 
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APPENDIX A: HIGH-ENERGY ASYMPTOTICS AS A SCATTERING FROM 

SHOCK- WAVE FIELD. 

The structure of the answer ( p5| ) for the high-energy scattering from external field can 
be made transparent if instead of rescaling of the incoming photon's momentum (|T^ one 
boosts instead the external field: 
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dx I dz6{z.)e'P^^T{j,{x + z)Uz)})a = I dx I dz6{zo)e'PA -{T{j^{x + z)Uz)})b 



(Al) 



where p^2^ = pf"^ + and the boosted field has the form: 



X 

Boi^Xo-! X±) = XAo{~r~ f x^X^ x±) 

A 

1 Xq 

A A 

5_l(Xo,X*,X_l) = A±(y,X=^A,X_L) (A2) 

where we used the notations Xo = x^p^i^, = x^p2^i- The field 

2 (f\\ 2 

A.(xo, X*, = AJ—XoPi H a;*p2 + (A3) 

So So 

is the original external field in the coordinates independent of A so we may assume that the 
scales of XojX* (and in the function (|A3|) are 0(1). First, it is easy to see that at large 
A the field B^{x) does not depend on Xo- Moreover, in the limit of very large A the field 
has a form of the shock wave. It is especially clear if one writes down the field strength 
tensor G^,^ for the boosted field. If we assume that the field strength F^^ for the external 
field vanishes at the infinity we get 

G^i{xo,x^,x_i_) = XF^i{^,x^X,x±) 6{x^)G±{x_i_) 

G.i(xo,x*,x_L) = jFj{^,x^X,x±) 

Xq 

A 

X 

Gik{xo, x^,x_i_) = -Fjfc(y , x*A, x_l) (A4) 

so the only component which survives the infinite boost is Fo± and it exists only within the 
thin "wall" near x^: = 0. In the rest of the space the field B^ is a pure gauge. Let us denote 
by fl the corresponding gauge matrix and by B^ the rotated gauge field which vanishes 
everywhere except the thin wall: 

B^ = limx^^-^GiSOAx*,x±) ^ 5{x,)-^G,{x±), B^ = B^ = (A5) 



Let us find the quark propagator in the B^ background (see Fig.|T^). We shall at first 
calculate the propagator in the external field B^ and after that make the gauge rotation 
back. 

We start the path-integral representation of a Green function in the external field: 
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FIG. 13. Quark propagator in the shock-wave field as a path integral. 



Pexp{igj^dt{B^^{x{t))x^{t) + \a^'' G%{x{t))} 

where a^y = ^i'y^'ju — lulfi)- First, it is easy to see that since in our external field (|A4|) 
the only nonzero components of the field tensor is G^^^ only the first two first term of the 
expansion of the exponent exp{Jdt^{(rG^)} in powers of (crG) survive. Indeed, o-^^G^^ = 
Yihl^^^i therefore {aG^Y ~ (l^27*)^ = since commutes with 7^. So, the phase 



factor for the motion of the particle in the external field ( |A^ ) has the form 



Let us consider the case x^: > 0,y^: < as shown in FigJT^. Since the external field exists 
only within the infinitely thin wall at x=k = we can replace the gauge factor along the actual 
path Xfj_{t) by the gauge factor along the straight-line path shown in Fig. ^ which intersects 
the plane = at the same point {zo, z±) at which the original path does. Since the shock- 
wave field outside the wall vanishes we may extend formally the limits of this segment to 
infinity and write the corresponding gauge factor as U^{z±) = [—ocpi + z±, oopi + z±] where 
label ^ reminds that we calculate this eikonal factor in the field B^. The error brought by 
replacement of the original path inside the wall by the segment of straight line parallel to 
Pi is Indeed, the time of the transition of the quark through the wall is proportional to 
the thickness of the wall which is ~ ^ which means that it can deviate in the perpendicular 
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directions iside the wall only to the distances ~ Thus, if the quark intersects this wall 
at some point {z^,z±) at the time r' the gauge factor (|A8|) reduces to: 



U'^iz^) + -^Ad,U''{z^) (A8) 



where the last term was obtained using the identity 



-— / dx^[oopf^ + x_L, —x^pf^ + x±]G^i{—x^pf'' + x±)[—x^p'i^ + x_L, -oopf^ + x±] (A9) 

So J So So So 



and the factor x^{t') in eq. ( [A^ ) comes from changing of variable of integration from t 
to x^^{t). Similarly, the phase factor for the term in the r.h.s. of eq. ([X^) which contains 
P^{x{t)) = j^'p^B^i^xir)) in front of the gauge factor eq. (|A6|) can be reduced to 

d 2 

-i)2T^[—X*pf+X^^-^ + XA] = -p2^{x,)[U{x±) - 1] (AlO) 
OXif Sq 

(The factor ~ (cG) is absent since it contains extra ^2 ih, ~ 0)- insert the 

expression for the phase factors (|A7| , |A10D into the path integral (|A6|) we obtain 

-i^,5(x.)[f/^(x^) - i]/o°°rfrAr-i/;Wr;px(t)e-^/>*^ - 

f S^dr rodr' Jdz6{z,)Af~' tiXlT^xit) /fc(r) 
e-/>4{f/n(^^) + _j_ /)t/^(^^)^JAr-i /:(Si'I^^Px(t)x.(r')e-^/>*^ (All) 



The additional Jacobian factor x^{t') in the numerator in the second term in r.h.s. of this 
eq. comes due to the fact that we must integrate over all r' from to r and therefore we 
insert 1 = /(ir'x*(r')5(a;=i,(r') — z^)) in the functional integral ( |A6|) . It is convenient to make 
a shift of time variable t' and to rewrite the eq. ( |A10|) in the following way: 
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- mx^)[U^{x^) - 1] CrfrAT-i St^:;Vx{t) Mr)e~'^r'^'' 
- h I^drS^dr' Jdz6iz^W-' 4)IT^^^(^) Mr)e-' St^Jy^^^t) 

{x,{t'W{z^) + I Af/^(z±>^2}e-^/.'''*^ (A12) 
Now, using the path-integral representations for bare propagators: 



J rx{T)=X . rT,^^2 1 

drAf-' I VxmT)e~^Io''^ = --- (A13) 

Jx(o)=y A-K\x-yy 



and 
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x{t)=x 



Vx(t)x^{r)e 



7r2(x — uY 



(A14) 



it is easy to see that the path-integral expression for the quark propagator in the shock-wave 
field ( |A13|) reduces to 



5(s,)[f/^-l](xx)+ \dzb{z,) 



27r2(x - z)- 



27r2(2-?/)4 

-i 



(A15) 



(in the region > 0,?/^, < 0). It can be demonstrated that the answer for the propagator 
in the region < 0,|/* > differs from eq. (|A15| ) by the substitution ^ U^'' . Also, 
the propagator outside the shock- wave wall (at , < or x,,, , > 0) coinside with bare 
propagator so the final answer for the quark Green function in the background can be 
written down as: 



(x 



y), 



2vr2(x-,)^+V^"^(^*) 
where we have used the formula 



^g^-M_|[f/n_i](,^)e(x.)e(-,.)- 



'. j dz5{zY} 



27r^{z - vY 



27r2(x - zY^^ 27r2(z - yY lix'^ix - y) 



;(e(x,)-e(y,)) 



(A16) 



(A17) 



to separate the bare propagator. In the momentum representation this answer ( |A16| ) takes 
the form: 



(A18) 



27rz(5(ap 



^ al(3kSQ - kl + ie sq 
0(-"°J(f/^^(p)-47r25(px))] 



which agrees with eq.(^) after integration over and rescaling ap = ja^ (here a° is the 
Sudakov component along vector p^). 
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Now, one easily obtains the quark propagator in the original field -B^ eq. ( [A2| ) by making 
back the gauge rotation of the answer ( |A16D with matrix It is convenient to represent 
the result in the following form: 



((. 



f 



(A19) 



+ I 
where 



27r2(a; - z^ 



2TT^{z-yy 



U {z±; X, y) = [x, z^] [z^, Zy] [zy, y] , z^ 



2 2 

{—ZoPi'^ H x^p2, z±), Zy = z^{x^ ^ y^) (A20) 

is a gauge factor for the contour made from segments of straight lines as shown in Fig. WA. 



(Since the field outside the shock-wave wall is a pure gauge the precise form of the contour 
does not matter as long as it starts at the point x, intersects the wall at the point z in the 
direction coUinear to p2 and ends at the point y). 

For the quark- ant iquark amplitude in the shock- wave field (see Fig 14) we get 




FIG. 14. Quark-antiquark propagation in the shock wave. 
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1 






1 


f 




{y 


f 



W{z^ 



where we can write down the gauge factor W{z±_; z'j_) = U{z±; x, y)W{z'j_; y, x) as a product 
of two infinite Wilson-hnes operators connected by gauge segments at ±00: 

W{zr,z'^) (A22) 
= hm [-upi + z±,upi + z±][upi + z±,upi + z\][upi + z\, -upi + z',][-upi + z\, -upi + zA 

The precise form of the connecting contour does not matter as long as it is outside the shock 
wave. We have chosen this contour in such a way that the gauge factor ( A23 ) is the same 
for the field and for the original field (see eq. ([A2|) . Now, substituting our result for 
quark- ant iquark propagation ( |A22| ) in the r.h.s of eq. ( |Al| ) one recovers after some algebra 
the eqs. (BSl BB) for the impact factor. 



APPENDIX B: ONE-LOOP EVOLUTION: WILSON LINES IN A SHOCK- WAVE 

BACKGROUND. 

Let us now find how the one- loop evolution of Wilson-line operators can be obtained using 
the shock-wave picture of high-energy scattering. To this end, consider the near-light-like 
operators U'^ and tl'^^ in the external field. Making the rescaling ( [A^ ) we obtain: 

{[OOPA + 2;_L, -OOPA + y±][-OOpA + X^, OOpA + = 

where the shock-wave field is given by formulas (|A2[ - |A^ ) We must find the derivative C,^ 
given by eq. (|69D. After rescaling according to eq. (^) one obtains: 

C^^{U{x^)U\y^))A 
2 

= ig— / udu{[oop'A + x^^up^^ + xa]F^o{up''^ + Xi_)[up^^ + -oop''^ + XA]U\y^))B - 

So J 
2 

i^f— / udu(p{xx)\-oop^^ + ?/_L, up^^ + yA\F^o{up^A + Vi^Xw^A + '^Pa^ + 2/±])b (B2) 

So 

Since the (=,,0) component of the field strength tensor vanishes for the shock- wave field ( |A4|) 
the only nonzero contribution comes from the diagrams with quantum gluons. In the lowest 
nontrivial order in ag there are three diargams shown in Fig. ^ 
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\ 



FIG. 15. Path integrals describing one-loop diagrams for Wilson-line operators in the shock-wave 
field background. 



Consider at first the diagram shown in Fig. |T^a (which corresponds to the case x,,, > 0, 
y^: < 0). The corresponding contribution to r.h.s. of eq. ( p^ ) is: 



{up'a +xa\up*{{Pa^ 



(0) 



1 



1 



P2p 



y±\ 
+ 



P2X 



ap 



(B3) 



As we discussed in Sect. 4 terms in parentheses proportional to Vo vanish after integration 
by parts (cf. eq. ([73|)). Further, it is easy to check that since the only nonzero component 
of field strength tensor for the shock wave is G^j_ the expression in braces in eq. ( |B3|) can be 
reduced to O^^ where the operator O^^, is given by eq. (|75|). Starting from this point it is 



convenient to perform the calculation in the background of the rotated field ([A5|) which 
is everywhere except the shock-wave wall. (We shall make the rotation back to field B in 
the final answer). Then the gauge factors [oo,M]t"[M, — oo] and [oo, f jt^fw, — oo] in eq. (p3|) 
reduce to t°[oo, — oo] ® t^[—oo, oo] (at > 0, < 0) and we obtain: 



- gH'^U'' ® t'f/^^ fdu f dv{u - v) (upj^ + x^\p,0^,\vp^^^ + 



ab 



(B4) 



where we have used the fact that the operator p^, commutes with O^. Let us now derive 
the formula for the (oo) component of the gluon propagator {x\0^\y) in the shock- wave 
background. The path- integral representation of {x\0^^\y) has the form: 



if dr 
Jo JO 



dr' i^x 



p2 



2p^ 



2p* 



Jo 



'2p/V' 
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rx{T)=X . fT , j.^ ri ri . rr 

drM-^ / Vx{t)e-'Jo'^'-{4: / dr' / dr'Te'^l' 



x(0)=y 



"'0 



p^ig J^„dtBj}ix{t))x^,it) r fl^iip^ig f^„dtBj}(x(t))x^it)gQn^^^^/r^^p^ig dtB^ (x{t))x ^,{t) 
<j 

JO a;*(r') 

As we discussed above, the transition through the shock wave occurs in a short time ~ ^ so 
the gluon has no time to deviate in the transverse directions and therefore the gauge factors 
in eq. ( P5| ) can be approximated by segments of Wilson hues. One obtains then (cf. eq. 



((. 



2 JO 

1 



^sl dr dr' dz6{z,)N'-' Vx{ty' 1''''- 



x{t')=z 



x{0)=y 



(B6) 



where [GG]*^ and [DG]^ are the notations for the gauge factors (^) calculated for the 
background field B^: 

[DG]^{x±) = J du[oopi + x±, upi + x±]D°'G^^{upi + x±)[upi + x±, -oopi + x±] 

[GG]^(2;_l) = Jdu J dvQ{u — v)[oopi + x±,upi + x±]Gl^{upi + x±) 

[upi + x±, vpi + x±]G^„{vpi + x±)[vpi + x±, -oopi + x±] (B7) 

As we noted in Sect. 4 the gauge factor —i[DG] + 2[GG] in braces in eq. ( P5|) is in fact the 
total derivative of U with respect to translations in the perpendicular directions so we get: 



O: 



2 JO 



/■cxj rr r /•x{t)=x ■ i ±^ I rx\T )=z . pi 

sl dT dr' dzS{z,)Af-' Vx{tyi^''''--—-dlU''{xA_)N~^ Vx{t)e-'Jr 



x{t')=z 



1 

xJr 



x{t')=z 



x{0)=y 



Using now the path-integral representation for bare propagator (|A13|) and the following 
formula 



drAf-^ 



Vx{t)^-e~'Jo = t )- — ^ 

x{o)=y x^{0) 167r^(x-?/), 



(B9) 



we finally obtain the (oo) component of the gluon propagator in the shock- wave background 
in the form: 



X 



O 



ln(x 



I y^^^(^*) 16^2^ [d^U''{z^)Q{x,)e{-y.)-d^U^''{z^)Q{-x,)e{y.)] 



(BIO) 



47r^(z — uY 



52 



where we have added the similar term corresponding to the case < 0,y^ > 0. We need 
also the ^ derivative of this propagator (see eq. ( P^ ) which is: 



ISr 



647r4 



dz- 



X - y)' 



1 



{z - yY 
(Bll) 



In the momentum representation this eq. takes the form: 



{k\0l\k-p 



ab 



-ISq 



2vr5K) — 



1 



(e(a°)(aif/^(p^)) - e(-a°)(9if/t^(p^)) 



(B12) 



aliPk- f3p)s-{k-p)l + te 



which agrees with eq. (^) after rescaling ap 
eq. ( p^) one recovers eq. (ffOD after some algebra 



jOp. Substituting now the eq. ( pil|) into 



y±) rU^{x±) ® t''W^{y±) + 

// ab 

y±]^U''{x^)t-®W''{y^)t'' 



(B13) 



Let us consider now the diagram shown in Fig. |l5|b. The calculation is very similar to 
the case of Fig. p!5| a diagram considered above so we shall only briefly outline the calculation. 
One starts with the corresponding contribution to r.h.s.of eq. (|B^ ) which has the form (cf. 

my- 



t^[vpf + x±, -oopf + x±] ® U^{y±) 

Crf + x,\up.{ipf, - - pi^^^^i 



P2X ^(3 
P-P2 



p-P2'V^gp^ + 2iG 



P-P2 P-P2 

(0) 



PV 



P-P2 " °^ 



As we demonstrated in Sect. 4 the terms in parentheses proportional to P, vanish and after 
that the operator in braces reduce to Ooo- Again, it is convenient to make a gauge trans- 
formation to the rotated field ( [A^ ) which is everywhere except the shock wave. Then the 
gauge factor \po^v\t°'\u,v\i^\v, — oo] in eq. (pi3|) reduces to ^^[oo, — oo]t* (at > 0,?/=k < 0) 
and we obtain: 

- g'tU^t'' ® [/^^ ^du^ dviu - v)lupf + xi_\p,0^,\vpf + x±)^^ (B15) 



Using the expression ( |B11|) for the gluon propagator in the shock-wave background after 
some algebra one obtains the answer (BHf): 



ap 



(B14) 
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x± 



ab 



(B16) 



The contribution of the diagram in Fig. [T^ differs from eq. (|B16|) only in change U ^ 
X ^ y (see eq. ( |A22|) ). Combining these expressions, one obtains the answer in the 
rotated field f|A5|) in the form: 



167r3 



dz I 



' [x - zjliy - z)\_ 
[{f/^(z^)};.Tr{f/^(x^)t/t^(z^)} - iV,{[/^(a;^)};.)f/t^(y^)f- 



1 



X — z) 
1 



(B17) 



Now we must perform the gauge rotation back to the "original" field B^. The answer 
is especially simple if we consider the evolution of the gauge- invariant operator such as 
Tr{[/(x_L) [x_L, ?/_L]_t/"''(?/x) [?/_L, + } where the Wilson lines are connected by gauge segments 
at the infinity, see eq. (P^D. We have then: 

C^(Tr{f/^(xx)[a;x,yJ„f/t^(y^)[l/^,a;^]+})^ 

= J dz±{Ti{U zj_]_U\zji)[z±,x±] + } Tr{U {z±)[z±, y±]-U\y±)[y±, zj_] + } 



iV,Tr{f/(xx)[x^,yJ_f/t(y^)[y^,xJ + })- 



ix± - y±y 



'x±_ - z±y{z± - y\ 



(B18) 



where we have replaced the end gauge factors like VL{oopi + x \_)VL^ {oC)pi + y _[_) and oopi + 
xj_)Vl} {—oopi + yA_) by segments of gauge line y±\+ and [a;_L, respectively. Since the 
background field is a pure gauge outside the shock wave the specific form of the contour 
in eq. (P18| ) does not matter as long as it has the same initial and final points. Finally, note 



that the gauge factors in r.h.s. of eq. (P18|) preserve their form after rescaling back to the 
field so we reproduce the eq. ( p.001 ). 

It is instructive also to see along which variable do the leading logarithmic integration 
actually goes. To this end we must find matrix element of the operator U{x_[_)U'^{y±) (see 



r.h.s. of eq. ( |B1|) ) in the shock-wave background. In the first order in one has (cf. eq. 



{U{x^)U\y^))Bn = -tgH^U''{x^)®t'U^''{y^) du di; (np^?) + 



O. 



(0) , 



ah 

(B19) 
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(we shall calculate only the contribution ~ which comes from the region > 0, < 
- the term ~ coming from x^, > 0, ?/* < is similar, cf. eq. ([B13D). Technically it is 



convenient to find at first the derivative of the integral of gluon propagator in r.h.s. of eq. 
( |B19| ) with respect to x±. Using formula ( |B11| ) for the gluon propagator (a;|C||/) we obtain: 



du / dv i^up^^^ + x± 



vp^A + yi) 



ab 



(B20) 



167r4 



dz I / — dv / dza 

u 







ab 



[u{uC,so - 2zo) - (x - - ie][v{vC,so + 2zo) - (x - z)\ 
The integration over Zo can be performed by taking the residue and the result is 

^2 >■ roodu. {x^-Z^UdlU{z^)] 



— I 



9 
167r^ 



dz^ / —dv — --"^'^"-L" y~^Aab 

u [{x — z)']_v + {y — z)\v — uv{u + v)C,sq + ie] 



(B21) 



This integral diverges logarithmically when m — O - in other words, when the emission of 
quantum gluon occurs in the vicinity of the shock wave, (note that if we had done integration 
by parts, the divergence would be at w — > so there is no asymmetry between u and v). The 
size of the shock wave z^ ~ m~^/X (where 1/m is the characteristic transverse size) serves 
as the lower cutoff for this integration and we obtain 



— I 



167r3 



In A dz 



1 da (x_L - z±)i[dlU{z±)] 



ab 



a [(x — z)']_a + {y — 



lnA((, 



X_L 



Pi,a2"^ 1 



yi. 



ab 

(B22) 



(recall that a = 1 — a). Thus, we have the contribution of the diagram in Fig. in leading 
logarithmic approximation in the following form: 



(f/(xx)f>Uy±))Bn = -(|^lnA)rt/^(xx)®f/t^(y±)((xx -,{dlU)-. 



y± 



ab 



(B23) 



which agrees with the first term in eq. ( |B13| ) (recall that ( 



1 X '9 ^ 



APPENDIX C: GLUON PROPAGATOR IN THE AXIAL GAUGE. 

Our aim here is to derive the expression for the gluon propagator in the external field in 
the axial gauge. The propagator of the "quantum" gauge field A"^ in the external "classical" 
field A'^'- in the axial gauge e^A^ = can be represented as the following functional integral: 

Gt{x,y) = limiV-i / nAAl\x)Al\y)e'I'^^^'''-^^^^^'^'''-^''^'-''^^^^^ (CI) 

where = — igAI^. Hereafter we shall omit the label "cl" from the external field. This 
propagator can be formally written down as 



55 



zGtix, y) = (x — ^ yV^ (C2) 

where = V'^g^'^ + 2igF^^ . It is easy to check that the operator in r.h.s. of eg. (|C2|) 
satisfies the recursion formula 

1 



ui 

which gives the propagator as an expansion in powers of the operator D\F^^ = —gipt°''yc,'ip- 
We shall see below that in the leading logarithmic approximation we need the terms not 
higher than the first nontrivial order in this operator. With this accuracy 

1 



Taking now w ^ we obtain the propagator in external field in axial gauge in the form: 

where the dots stand for the terms of second (and higher) order in D'^F\p. It can be demon- 
strated that for our purposes a first few terms of the expansion of operators ^ in powers of 
F^ri are enough, namely 

^G%{x,y) 
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